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Final Report to Innovative Research Program 
Precision Optical INTerferometry in Space (POINTS) 


This is the final report, and incorporates the several semiannual reports listed on the 
cover, of the work done under NASA Innovative Research Program Grant NAGW-1647 for 
Precision Optical INTerferometry in Space (POINTS). During the tenure of the grant, the 
POINTS concept was developed and supporting technologies were investigated. One 
technology of particular importance to POINTS is laser gauging. Under the grant, we 
developed a new class of precision laser gauges. In our laboratory, we demonstrated a laser 
gauge at the picometer level (10'*" meter) required for the POINTS spacebome instrument. 
This is a three order advance over the best commercial gauges available. A patent is pending, 
and industrial applications are being sought. 

Appended are two recent invited papers on POINTS. The paper to appear in the 
special issue of Remote Sensing Reviews provides some historical perspective on the POINTS 
project, underscores the innovations, and indicates the work by students on the project. The 
paper for the SPIE Conference 1947, Spacebome Interferometry , contains the most up-to-date 
description of the instrument. Between them, these papers summarize the status of the 
development at the end of the subject grant. (The SPIE Conference included five other 
papers on POINTS as indicated in the copy of the Conference agenda and the POINTS 
bibliography that follow. Copies of the papers in the bibliography are available on request.) 

POINTS is now a collaborative effort between SAO and JPL. The development of 
POINTS at SAO is now supported by the PIDDP program of Code SL and by the 
Smithsonian Institution. We plan to propose to further develop POINTS in response to the 
NRA that is expected to be issued this summer for the TOPS-1 mission of Code SL. 
According to the agreed division of labor, JPL will buy or build the spacecraft and do the 
instrument-spacecraft integration; SAO will develop the instrument with the help of industry. 
After an evaluation of potential partners, we have selected Itek of Lexington, Mass, as our 
industrial partner for the instrument. 

Also appended is a recent paper on a small astrometric interferometer named for 
Simon Newcomb. This is a spinoff of the POINTS work and is now the subject of a 
collaboration between SAO and the Naval Research Laboratory. Newcomb is intended for 
the Navy’s Space Test Program, which is a fast-o*ack low-cost progiam foi technology 
development and evaluation. Under the established division of labor, SAO will be 
responsible for the instrument. Itek is designated as our industrial partner for the Newcomb 
instrument. NRL will build (or possibly buy) the spacecraft and do the integration. The 
Space Test Program will provide the launch. 
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POINTS: the first small step 


R.D. Reasenberg. R.W. Babcock, M.C. Noecker, and J.D. Phillips 

Smithsonian Astrophysical Observatory 
Harvard-Smithsonian Center for Astrophysics 


ABSTRACT 

POINTS, an astrometric optical interferometer with a nominal measurement accuracy of 5 microarcseconds for 

OSSA 8 IwlirjZrf 1^ SeP T!? by f^° Ut 90 dCg ’ ‘ S present,y under consideration by two divisions of NASA- 
^ ^ powerful new multi-disciplinary tool for astronomical research. If chosen as the TOPS-1 (Toward 

“ Systems) instrument by the Solar-System Exploration Division, it will perform a definitive search for 

extra-solar planetary systems, either finding and characterizing a large number of them or showing that they are far less 

^T^“ 1,eved - ? Ch0SC " aS r AIM (AStr0metric I^rferometry Mission) by the 

INTS wiU open new areas of astrophysical research and change the nature of the questions being asked in some old 
s. n either rase, it will be the first of a new class of powerful instruments in space and will prove the technology for 
the larger members of that class to follow. Based on a preliminary indication of the observational needs ofTe ^o 
missions, we find that a single POINTS mission will meet the science objectives of both TOPS-1 and AIM. 

in „ , ThC de,cc,s a dispersed fringe (channelled spectrum) and therefore can tolerate large pointing eirors 

°|*!f ,0n K . the dlf fi cul1 Problem of measuring the angular separation of widely spaced star pairs is reduced toLo less 

it^mlti Pr< ff e , mS f L° f m ! QSUnn u g ,he 311816 belween the lwo interferometers and that of measuring interferometrically 
the small offset of each star from the corresponding interferometer axis. The question of systematic error is the central 
Uieme of the instrument architecture and the data-analysis methods. Stable materials, precise thermal control and 
continuous precise metrology are fundamental to the design of the instrument. A preliminary version of the required 
»ZZ;\ S" mC ! r h 0 2y haS I™ dcmonstra,ed in the laboratory. Post-measurement deletion and correct!™ of time- 

tm^fnnc b w C CmC " ,S in da,a a,1alysis - ln thal post-measurement analysis, individual star-pair 

separations are combined to determine both the relative positions of all observed stars and several instrument parameters 

time-dependent measurement bias. The resulting stellar separation estimates are both global and bias- 
free at the level of the uncertainty in the reduced (i.e., combined and analyzed) measurements. 

1. INTRODUCTION 

New technologies have made possible construction of ground-based and space-based optical interferometers that 

working C 0b , SerVa,, °" S ° f hi8h ***** value - Ear» y versions of the ground-based instrumms Tnow 
^'" gOF n6 f completion. However, ground-based interferometers and telescopes jointly suffer two classes of 

bS f m ?° n : “* ma<JC thr ° Ugh 3 ,Urbulenl almos P here wi « b limited windows of visibility; and (2) The 
platform (Earth) rotates slowly, necessitating compensatory pointing, restricting the portion of the sky that ran be 

° bS !^ d 31 0n ? Um6 ’ and seIIing ^ ,ime between observations of points well separated on the sky. These restrictions 
die Whe " ‘ he 0bservall0ns ^e made from space, where several optical interferometers ran be expected during 

the 2lst century. We expect that most of the starlight interferometers in space will be imaging devices ona grand scale 
with corresponding pnee tags. Because they represent a new and complex technology, these "greater observatory" class’ 
mshuments are not likely to be flown without the benefit of a smaller, less expensive^ecuraoHnstrament^as 

POINTS, a space-based astrometric optical interferometer with a nominal measurement accuracy of 5 
microarcseconds (pas), will be a powerful new multi-disciplinary tool for astronomical research. 1 ' 2 The instrument 
includes two separate stellar interferometers thal have their principal optical axes (nominal target directions) separated by 
9 - an adjustable angle of about 90 deg. (See Fig. 1.) Since the targets are within a few arcsec of thei resSvT * 

axes ; 0ff ' axi l dis,0rtl0ns ,he a,lendant biases are minimized. In each interferometer, the starlight 
Iraving the beamsplitter is disperaed to form a "channelled spectrum." which gives the instrument an enhanced tolerance 

fj™ ofmr^l 6Spe i ,ally hel P ful dunn = ***« acquisition. A single measurement determines the angular separation of 
Thp r n • g f arS ’ ^ acb mterferometer has a baseline 2 m long, and two afocal telescopes with 10:1 compression 
The two interferometric subapertures of each starlight interferometer are 25 cm in diameter. Our studies of detectors and 
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350 - £* - * if - -»“ - ™ 8 10 - 
r .^^sssssE^ur^^ssi^tsJS 

^^ b r«fc™" S Sen IL .he 

Of finding a suitable reference. If the relative ^ tcd (0 ^in about 80 stars as bright as (visual) 
has an area of 2160 square e F c< ^ c , Q Thus the observation time is not dominated by the low photon 

mag 5; 1200 stars, mag 7.5; and ’ is 8 a care ^i y ’ studied target star, chosen from the reference grid described 

rsSS*" mTmotan, of .he reference erne are well modelled and do no, connp, the measurement of oO>er 
stars. 

The instrument uses two kinds of metrology systems. The first is the full aperture metrology (FAM), which 

S5S^=S3SSr- 

interferometers. 

In FAM modulated laser light is injected backward into the starlight optical path and fully illuminates the optical 
starlight interferometers. 




J3H=Sr25S.=«; 

Schumaker et a/. 10 


Reasenberg 2 




1947-02 


3 ”‘,° f f ■ <” lm "°-P°™ ^ Itoges tha, measure 1,0* six 
measured. The small diSce ta»Sl ,nleres ‘i ““ a " 8 " pseudobaselines is 

from the nominal optical axes of the starlight interfemmete fJ" f econ ^ order on the displacements of the targets 
pseudobaselines. Thus, the measurement of one an<>le vieldk ?h . ,ween ,he baselines ar) d the corresponding 

bias are discussed in Section 2 4 2 ^ ^ 71,6 parame,ere of ^ 

in the astrometric data, as discussed in Section 4. PP ^ X ' C b,aS * S esl,ma,ed from 1,16 closur e information 

Mission 


INSTRUMENT DESCRIPTION AND CONCEPT 


2.1 Principles of operation 


angle, and a SC Wot* * “ arly ^. 8h ' 

sSStoSS S l,rin8 f ndq,end ' n ' ,y s ‘ and 5: ' lhe <«*£*£ 4? 

<*- P-. Which uses 

toward £ S^ra^S^tblvC 23“ 7TT ^ « * T**- -*« - <tocto 
onto a pair of linear arrays of detectors. When the star is on axis (8 - "" l ' gl “ K dlspetsed ” d totused 

intensity at the two beamsplitter exit pons When the star is off Jkffi !' m 8 " a * any S "' e " wavelen glh has equal 
wavelength a, one pon is complement by dS£. imerf^ „£ auhe "T*™' «* a »«n 

r;7n,“:rvC s ^ 

^ whi,e S'" frinSE d£KC,ion would ^ a " 

malting the initial acquiS/X 2? “ ^ “ 2 

pointing system and removes the need for ,h°L, ' ' he USe ° f ,he dls P ersed fnn 8 e greatly simplifies the 

add mis Ld complexily > and red^MUtv tL ° conTI/^ d “™ 8 a " »fservation. whtch would 

additional informaL availahk in L ctaCL^mm ?? ° f d ' SPmed frin8e a PP roach is «* 'toe -s 
Closely spaced and '° ^ “ 8ets "> a ' « 

technique exclusively ias preposed b, S. J Endal^ ™ ^ ^ M mS ' rumem “ s 

2.2 Detectors and instrument pointing 

• .w. ThCre ^ threC quesIlons concerning detectors and pointing: (1) How far from the target direction cem »h 
instrument point and still delect the target? (2) What is the effect on the inform? , V g direct,on can die 
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The information content of the data is monotonic with V. both of which are zero at 6=28*. where there is one fringe per 
pixel, nie information content falls rapidly with decreasing V near V = 1, so the viability must be kept high. A pnsm. 
like a diffraction grating (but in the opposite sense), provides a nonlinear dispersion so that, if the detector pixels are o 
untform' wkhh? the number of pixels £r cycle of channelled spectrum will vary with color. One solution is to combine a 
grating and a prism, as is being done for the IOTA ground-based interferometer. 6 An alternative approach, which is 
adopted for POINTS, is to use only a prism, allowing the pixel width to vary with the color and so avoidrng the hght loss 
of the grating. Below, we consider in turn the detectors, the spectrometer, limiting magnitude, and the fine pointing and 

isolation system. 

As discussed below, we have found that the fringes can be held sufficiently stable on the detectors that an 
integrating detector such as a CCD can be used with little loss of fringe visibility due to motion during the mtegrauon 
period. CCDs offer high quantum efficiency, and there is a wide base of experience of their use in space as well as 
continuing development. With a stabilized fringe integrated on the CCD detector, there is an automatic solution to the 
problem of data compression, which otherwise would require considerable on-board computation. As presently 
envisioned, the CCD P arrays will have 650 cells to cover an optical bandpass from 0.25 to 0.9 pm, and the nominal mirror 

coating will be aluminum. 

The spectrometers, one for each exit beam, must disperse the starlight to produce the channelled spectrum while 
maintaining high fringe visibility. Three primary effects govern resolution and thus fringe visibility: 
aberration, and detector pixel size. In the dispersion direction, the diffraction-limited spot must be £ 1/30 the spacing o 
the dispersed fringes in order to maintain high fringe visibility. Additionally, providing spatial resoluuon in the cross- 
dispersion direction might allow simultaneous acquisition of data on multiple targets in crowded fields. To improve the 
performance in crowded fields, a slit could be provided to eliminate light from stars too far off-axis to produce 
measurable fringes. This would require additional reflections and so reduce the single-star throughput. It has not yet 
been decided whether to include a slit. 

So far, we have considered slitless designs in which the collimated beam from the starlight beamsplitter passes 
through the prism and is focussed onto the detector by a single off-axis mirror. These designs are characterized by t ec 
parameters: focal length, prism deviation range, and pixel size. A range of designs meet the resolution requirement, 
we are currently determining what design best meets a range of secondary requirements including Nyquist angle 5*. 

cross-dispersion resolution, dark current, sensitivity to detector position, the option to include a polarizing spectrometer, 
and suppression of scattered metrology light. 

Limiting magnitude depends on dark current. We believe that present CCD's impose a limit somewhere between 
magnitudes 17 and 20. The two important factors for setting the limiting magnitude are detector temperature, current y 
exacted to be around -80° C. and the area of detector in which the stellar spectrum appears. 

In the present design, the instrument is connected to the spacecraft scan platform by means of a soft single-axis 
pivot with a free swing range of several tens of arcsecs. A flexural pivot that would serve well as thebas.s for die 
Isolator has been made by the Perkin-Elmer division that is now Hughes-Danbury Optical Systems (HDOS). (The pivot 
was flown as part of the Apollo Telescope Mount on Skylab, 1973-1974. 7 - 8 ) There are several suitable magnetic 
actuators to provide fine pointing correction to the isolated instrument. The instrument mounted on the pivot would have 
a torsional resonance between 1/20 and 1/5 Hz. b Given such a soft mount, there would be considerable filtering at the 
delector sampling rate of 100 Hz. applicable to bright target stars. 


• Consider an observation in which 10 7 photons are collected in 100 batches by one interferometer Without 
compression two numbers (epoch and detector number, about ten digits total) must be transmitted to the ground for each 
photon. Without stabilization and on-chip integration, the fringe must be rotated and integrated in “ ftware ; r ^"8 
several floating point operations per photon in addition to the real-time filter. With stabilization and on-chip integration, 
there is little computation required for fringe tracking, and the telecommunication load without explicit compression is 
2 x 10 4 numbers of five digits. It appears that a compression of at least 10 fold will be possible. 

b There would need to be cables, etc. crossing the pivot, and these would add to the stiffness by an as yet 
undetermined amount; this is the principal uncertainty contributing to the stated possible range of the frequency. 
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•r» fi h #h h K schem f [ equires *a« al leasi one target be bright, say mag 10. Initially, a star tracker would be used 

to find the bnght target and the approximate rotation rate of the instrument. This information would be used to reduce 
the offset and rotation rate to below the threshold for detecting starlight fringes. Once the fringes were found, they would 
serve as the reference for the fine pointing. On a mag 10 star, with a sample rate of 100 Hz, the samples would have a 
statistical uncertainty of »300 pas. However, the use of one of the interferometers as the angular reference requires real- 
time data analysts and the corresponding on-board computational capability. The 100 Hz sample rate probably limits to 
*IU Hz the unit-gain frequency of the control loop for the orientation of the instrument. 

23 Internal measurements and systematic error 


The ^ 0nt ^ 1 0f s y $tematic error is central to achieving the stated instrument performance and the mission science 
H S .' add T SS thls P rob ' em al three ,evels: (1) stable materials and thermal control; (2) real-time metrology; and 
(3) the detection and correction of systematic error in conjunction with the global data analysis. The last of these is 
addressed in Section 3.4 The best materials fail by orders of magnitude to provide the long-term dimensional stability 
required to maintain each interferometer’s OPD and other critical dimensions, in the absence of other means of control. 
Stable materials for the structural elements of the instrument serve to limit the dimensional changes that the metrology 
system must determine, and to control errors that are second order in component displacements. In a few places we are 
orced to rely on material stability (over short times). The instrument is designed so that such metering elements* are 
small* and can be thermally isolaied and regulated. 

For a 2-m baseline, the nominal 5-pas uncertainty corresponds to a displacement of one end of the interferometer 
toward the source by 0.5A = 50 pm (picometer = 10 1 - m). Since similar displacements of internal optical elements are 
afeo important, the instrument requires real-time metrology of the entire starlight optical path at the few picometer level. 
This metrology does not pose an overwhelming problem for two reasons. First, the precision is needed only for 

averaging times between 1 and 100 minutes. Second, a slowly changing bias in the measurement is acceptable as 
discussed in Section 3.4. H ’ 


The instrument relies on two kinds of laser-driven optical interferometers to determine changes in critical 
dimensions. New concepts and technology have been required for these interferometers to obtain picometer accuracy 
The aser gauges per se a new application for a HOE, and the laser-gauge endpoint assemblies were invented during the 
development of POINTS. The laser gauges can be grouped according to whether they measure distances internal to a 
starlight interferometer or distances between starlight interferometers. Below, we consider the former. Jn the next 

section, we consider the latter and some related questions of instrument geometry (the relation between the baselines and 
the pseudobaselines) and measurement bias. 


The high-precision star position measurement is made with respect to the optical axis of the interferometer bv 
measuring the difference in the optical paths from the target to the beamsplitter via the two sides of the interferometer. 

In turn, the position of this axis is determined (defined) by the positions of the optical elements used to transfer the 
7 he J m 1 etro i l i °gy s .y sIem mus ‘ determine, to about 10 pm overall accuracy, the average change in the starlight 
OPD s induced by all motions and distortions of all optical elements. Our approach is to use FAM, a novel technique 
discussed below. FAM provides three significant advantages over conventional approaches. 

(a) FAM removes complexity. The usual metrology systems use a large number of laser gauges to determine the 
locations of the elements individually. From these measurements, the optical path through the system is computed FAM 
directly measures the optical path through the system. y 


(b) FAM measures the correct quantity. Because the metrology signal fully illuminates the surface of each 

optical element that determines the starlight OPD at the beamsplitter, the metrology OPD represents accurately the 
average starlight OPD through the system. " 3 

(c) FAM provides the basis for an operational definition of the direction of the interferometer baseline It 
measures the OPD with respect to a pair of fiducial points that define the pseudobaseline, which is located in front of 
each interferometer and which is held at a small fixed angle (nominally zero) to the real interferometer baselines. These 
fiducial points are used to determine cp. the angle between the two interferometers’ optical axes. (The small bias 

resulting from the angle between the pseudobaseline and the real baseline is determined routinely as part of the data 
analysis.) v 
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Figure 2 illustrates the preferred version of the technique. The key element is the set of primary mirrorswh.ch 
have shallow (phase contrast) zone plates on their surfaces: alternate zones are depressed about 100A. In this HOE^ the 
zones are approximately in the form of Newton’s rings; each zone has about the same total area. For the flight hardware, 
we would manufacture the HOE by a method that would provide the required spatial dependence of diffracUve efficiency. 
The needed binary pattern would be found numerically and used to make a mask such as is used in the manufacture of 
integrated circuits. The mask could be written on a master plate by an "e-beam machine or a laser lathe. Each minor 
would be given a thin (ca. 100A) coating of metal and contact printed using photoresist and metal etching. Finally, the 
optic would be given a reflective coating. This technique has been used for many years at Hughes-Danbury Optical 
Systems, Inc. (HDOS) in support of DoD projects. 

It is important that the thickness of the metallic layer on the optic match the desired radial profile to perhaps 
Vi% as the diffractive efficiency is proportional to the square of zone thickness in the intended regime. The radial 
variation is required to compensate for the intensity variation of the metrology beam, which is roughly Gaussian To 
achieve the required accuracy, it may suffice to deposit material under carefully controlled condiuons, with masks to 
achieve the desired radial variation. A densitometer could be used to measure the thickness variations, asthe coating wi 
have an optical depth of order 1. It may also be possible to use Plasma-Assisted Chemical Etching, an HDOS process, to 
correct the profile after deposition. 

In the traditional FAM system, modulated laser light is injected backward into the starlight optical path and fully 
illuminates the optical elements that transfer the starlight. The hologram diffracts about 2% of the FAM light, and t e 
samples of FAM light from the two sides of the interferometer are brought together at the metrology beamsplitter to 
produce an error signal that drives a null servo. In the process, the fiducial blocks are sun-eyed into position. Each 
fiducial block contains four incomplete hollow comercube retrorefiectors. constructed such that then apices coincide at 
the fiducial point to within a few microns. Knowledge of the six distances defined by these four fiducial points 
determines the angle between the pseudobaselines of the two starlight interferometers. There is a substanual advantage to 
making the FAM laser light travel in the same direction as the starlight. We are investigating the principal disadvantage 
to this variation, that laser light will be scattered into the starlight detectors. 

The FAM technique requires two principal servos. The laser beams diffracted from the HOE s in the two arms 
of the interferometer are recombined at a metrology beam splitter; this is a Mach-Zehnder interferometer configuration, 
which is sinusoidally sensitive to the length difference from beam splitter to beam splitter via the two paths A servo 
regulates this path difference to a constant value by moving the starlight beamsplitter assembly. A second laser beam 
injected into the metrology beam splitter parallel to the FAM beams similarly measures the difference in distances from 
that beam splitter to each of two comer cube retrorefiectors situated in the fiducial blocks; this auxiliary servo regulates 
that path difference to a constant value by moving the auxiliary beamsplitter assembly. With both servos working there 
is a small (under 1 mm) and constant difference between the distances from the starlight beamsplitter to the fiduci 
points, which are the apices of the retrorefiectors inside the fiducial blocks. These servos can have small bandwidths^ 
because the time scale for distortion is long compared to a second, and because we require that vibrations of the OPD be 
small enough that they need not be tracked. Since it is desirable to limit the contamination of the starlight by the laser 
signal (e.g.. via scattering from the surfaces of the optical elements), a minimum of laser light is used m the FAM servo. 

We have shown that a change in the distance between the telescope primary and the atheimal collimating lens 
has an effect on the FAM OPD that is about 0.7% of the change of distance. For this reason, we introduce a third servo 
related to FAM that moves one (or both) of the fiducial blocks. (Note that if both fiducial blocks are moved in equal and 
opposite directions, the two overlapping illuminated regions on the auxiliary beamsplitter move in the same direction.) 

The sensor for this focus servo is a laser gauge that measures the difference of the distances from the auxiliary 
beamsplitter to the telescope primary miiTors. With the auxiliary servo working, the focus servo measures the dirierenu^ 
distance from the fiducial points to the primary mirrors. Note that the fiducial point is severe cm from the athermal lens. 
This distance, which must be stable to be about 0.1 nm, is metered by the body of the fiducial block, which is quite 
stable, as discussed below. 
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Figure 2. Inierferometer optical paths, (left) Starlight and auxiliary null interferometer. For the latter a laser signal 
passes through the spatial filter A and is divided by the metrology beamsplitter B to form beams that enter the fiducial 
blocks C and C . Within each fiducial block, the beam is deflected by the 45° minor, retroreflected by the hollow 
comer cube, and returned toward the metrology beamsplitter by the 45° minor. The returning signals are combined at the 
metrology beamsplitter and fall on detectors D and D'. The metrology beamsplitter may thus be refereed to the 
comercube apices, (nght) FAM interferometer. A laser signal passes through the spatial filter E and injection 
beamsplitter F and is divided by the starlight beamsplitter G. The separated beams are reflected from the tertiary minors 
H and H and secondaiy minors J and J ' to fully illuminate the primary minors K and K ' . The zone plates on the 
pnmanes diffract a small proportion of the metrology light to focal points in the holes in the secondaries, and the 
athermal lenses collimate the diffracted light. Finally, the signals are reflected from the 45° minors, are recombined at 
the metrology beamsplitter, and fall on detectors L and L 1 . The starlight beamsplitter may thus be refereed to the 
mettology beamsplitter via the starlight optics, and thence to the comercube apices. Note added in proof: In the present 
esign, the FAM light is injected at the auxiliary beamsplitter and travels in the same direction as the starlight See the 
paper in these proceedings by Noecker et ai. A 


Reasenberg 7 


1947-02 


The fiducial blocks (Fig. 3) and the metrology beamsplitter assemblies pose the only critical materials problem 
identified. To minimize changes of the size or shape of either of these, which would cause a corresponding changes in 
the bias of the metrology system, they will be made of a stable material and kept in a thermally stable environment. For 
the fiducial blocks, we plan optically contacted Premium ULE®; the beamsplitter material is still under study. The 
effective separation of metrology links within the fiducial blocks would be under 10 pm. and the temperature would be 
heVd suble ^ deg C on a time scale of a few hours. The effects of tempera.uie fluctuations on a longer ume scale 
will ^ routinely removed during the analysis on the ground of the astrometric data. TJe thermal gradient : wd be reduced 
by putting a pair of nested aluminum cylinders with 1.5 mm walls around the fiducial block. We plan to ach^e the 
required thermal stability with active control of the outer skin of the fiducial block enclosure to eg ( g 
deg C does not seem difficult), and with polished facing surfaces on the outside of the inner cylinder and ^ mside 
outer cylinder. The poor radiative heat exchange between the cylinders, combined with the hea “capacity _ of ^th > ™ucial 
block and inner cylinder, will yield a thermal time constant of about one day. Using Premium ULE , which has a 
coefficient of thermal expansion |a| < 5 x lO ’/deg C in a range from 5 to 35 deg C. the resulting distance error would 
be < 10' 3 pm. More important would be the effects of thermal gradients of up to 10 deg C over the - 5 cm separation 
of metrology beams where they reflect from the fiducial block mirrored surfaces. Such gradients would produce a 
measurement bias under 1 pm, again within the required limit. 

In the manufacture of the fiducial blocks, the critical concern is for the alignment of the "retrostnps," which are 
sections of hollow comer-cube retroreflectors. (See Figure 3.) Each fiducial block has four retrostnps. each of which 
must be assembled onto the main body from three small pieces of glass, using optical contacting, glue a few pm thick, or 
wettino with a thin layer of indium. Each retro has a symmetry axis that must point in approximately the correct 
direction and all four apices must be collocated to within a few microns. We now envision a manufactunng jig that 
manipulates and holds the three small pieces of glass in position on the main body. Dunng assembly the jig wou d 
permit the unit to be rotated on precision bearings, perhaps air bearings. The alignment procedure calls for the retrostnp 
under assembly to form one part of an interferometer using laser light. Under rotation, the changing fringe pattern wou 
indicate misalignments of the three small pieces of glass. Algorithms have been developed for converting the 
interferometer fringe pattern into corrections to be made to the positions of the three small pieces of gtoss. T e 
interferometric measurement is not expected to limit the accuracy of the alignment of the retrostnps_ The precision 
bearings, which should have sub-micron runout, and the devices used to position the small pieces of glass will prob y 

set the limit to the alignment. 


2.4 Determination of <j> 

The angle <p between the pseudobaselines of the two interferometers (/.<-.. the instrument articulation) is 
determined by the measurements of the six distances among four fiducial points in the system. 

costp = (d 2 + d 3 2 -d, 2 -d 4 2 ) /2b m b f 


Here d are the four distances between the fiducial points of one interferometer and the fiducial points of the other and 
b . b f are the baseline distances for each interferometer, /.e„ each is the separation of the fiducial points of one 
interferometer. Each of the six distances is measured by a laser gauge. (The laser gauges per se are discusse y 
Noecker el al., A who note that the techniques we have developed easily extend to an absolute distance gauge.) We have 
investigated the problems associated with the laser gauges reading the distance correctly except for a fixed bias, as wou d 
be the case if an incremental gauge were used and the zero point determined crudely by alternative means (e g., by an 
encoder on the articulation axis or from the known positions of bright stars observed by startrackers mounted on the _ 
interferometer optical benches.) Under these conditions, there would be an additional five parameters of the system that 
would need to be determined from the astrometric data each time the laser gauges were turned on. ^minal^ tas * ou ‘ d 
be at the start of the mission only. Our covariance studies show that estimating these parameters poses no problem even 
if it were decided to estimate them ob initio as often as every day. 

The metrology system measures the angle between the pseudobaselines, which are defined by the pairs of 
fiducial points in the two starlight interferometers. The two (nearly) perpendicular pseudobaselines define the 
"interferometer plane." For simple analyses, we assume that each starlight interferometer lies in the interferometer plane. 
However we have considered the effects of the three small rotations of each of the starlight interferometers. (1) arou 
the perpendicular to the interferometer plane; (2) around the pseudobaseline; and (3) around the starlight mterferometer 
nominal optical axis. The first type of rotation yields a bias proportional to the difference of the rotation angles. This 
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Figure 3. A fiducial block. There are four "retro-strips." each a slice of a truncated comercube retrorefiector. and each to 
serve as the end point for a laser gauge. The apices of the comercubes coincide at the fiducial point in this computer- 
generated engineering sketch. v 


bias is directly observed and held stable by the laser gauges and then estimated and removed in the data analysis. The 
other four rotations (two for each starlight interferometer) do not contribute to the measured angle in first order, but do 
contribute in second order. Thus the latter two types of rotations need to be aligned to somewhat better than an arcsec. - 
or further discussion, see Appendix. In a related study, we are considering how the above rotations and skews could 
come about from misalignments of the optical elements of the starlight interferometers. 

Although the above described metrology system is capable of providing the required precision, it contains many 
finite-size optical components, each of which will introduce a bias into the measurement of the angle cp. This bias is 
likely to be time dependent at the microarcsec level. It is essential that we be able to determine and correct for the 
instrument bias, preferably without the introduction of additional hardware. Because POINTS, like HIPPARCOS does 
global astrometry, that bias determination and correction naturally occur when the observations are combined in a 
least-squares estimate of the individual stellar coordinates (including proper motion and parallax), instrument model 
parameters and the expected biases. In particular, our covariance studies have shown that even without the introduction 
ot a special observing sequence, it is possible to estimate simultaneously the stellar coordinates and tens of instrument 
bias parameters per day. The stellar coordinate estimate uncertainties increase by a factor of 1+0.003 N where N is the 
number of bias parameters to be estimated per day, and M = 5. (The redundancy factor, M, is defined in Section 4 1 ) 
Thus, metrology biases and related errors can be allowed to change on a time scale of hours without significantly 
7 ,™ 8 * he of lhe instrument. This result allows us to relax the required stability at time-scales longer 

than 100 minutes of both the laser used in the metrology and the rest of the instrument. This subject is also discussed in 
Section 4. 


3. SPACECRAFT DESIGN 

The POINTS spacecraft will contain the usual set of systems to support the instrument: telecommunications, 
altitude control, power, and computation. Pointing interferometers at each of two stars requires control of rotation about 
four axes. One axis is provided by relative articulation of the interferometer axes. The instrument will be articulated 
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with respect to the solar shield with a two-axis drive. That 
articulation will be constrained by the requirement that the 
Sun shield keep the instrument housing (and preferably also 
the bus) in shadow at all times. The fourth axis is provided 
by rotating the spacecraft-plus-instrument about a line to the 
Sun. The spacecraft bus may require an isolator to keep jitter 
from entering the instrument from the spacecraft. For a more 
detailed description of the spacecraft see Schumaker et al. 

Within the instrument, there is an articulation 
mechanism that sets <p. This mechanism must permit about 
0.1 radian of motion and should have a resolution of order 0.1 
arcsec: 2 x 10 s steps. Of course, during observations, the 
combination of the articulation pivot and drive must be stiff in 
all six degrees, especially in the rotation that changes <p. 

The instrument is kept behind a solar shield, which is 
important in permitting large coverage of the sky without 
exposing the instrument interior to the direct light of the Sun. 

The solar cell array is mounted on portions of that shield. 

During launch, portions of the shield are folded along the 
spacecraft to fit inside the launch shroud. They provide 
limited power during the 17 hour cruise from parking orbit to 
the 100,000 km altitude at which the orbit is circularized. 

The question of sky coverage is normally considered in the context of an instrument with a single pointing 
direction and exclusion anales for the Sun. Earth, and Moon. For POINTS, the problem is more complicated because the 
instrument needs to look simultaneously in two directions about 90 deg apart. Thus it lakes four angles to describe the 
phase space of the sky coverage. For the purpose of calculating sky coverage, we constrain the target separation to be 90 
deg, reducing the number of independent angles to three. We may neglect both the Earth and Moon exclusions since we 
are willing to wait a few hours for any particular measurement. (The nominal orbit is circular at 100,000 km, which has 
a four day period.) Thus we are left considering only the spacecraft and the Sun. The line between them is an axis of 
symmetry. Rotation about that axis can be neglected in discussion and analysis of sky coverage; there are only two 
nontrivial angles. 

In our sky-coverage study, the first parameter was cos(y), where y is the angle between the anti-Sun direction 
and the first star, as seen from the spacecraft. The second parameter was w, the Sun-starl-star2 angle on the celestial 
sphere centered on the spacecraft, i.e.. the angle between the planes defined by: (a) spacecraft. Sun, and starl; and (b) 
spacecraft, starl, and star2. The study maps out the portion of phase space that is accessible. With this choice of 
parameters, the area in the "cos(y) - co plane'' is proportional to the volume of phase space. With the nominal 
configuration. 83% of phase space is accessible at any one time, and stars in 99.1% of the sky can be observed, but some 
with a limited range of co. (Stars in 98.5% of the sky can be observed with a cumulative range of co of at least 180 deg.) 
In all cases, the stars or star pairs that cannot be seen at a particular time can be seen a few weeks earlier or later. 

The original work has been extended to address the relationship between the available volume of phase space 
and the size of the solar shield. As shown in Fig. 4, increasing the size of the shield increases the sky coverage for 
shields similar in size to the nominal, which has a 30 m 2 area. However, saturation of sky coverage is found for shields 
with diameters larger than about 7m. 
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Figure 4. Fraction of star pairs with 90° separation 
visible to POINTS at any epoch. 


4. MISSION OPERATIONS 

Here we consider five aspects of mission operations: normal operation of the instrument, the effect of a faint 
companion, data analysis, the search for other planetary systems, and in-orbit testing. We do not discuss detailed target 
selection since there is an overabundance of interesting targets and the actual selection will depend strongly on which 
division of NASA/OSS sponsors the mission. 
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4.1 Normal instrument operations 

Our understanding of the characteristics of the 
instrument s operation is based on a series of covariance 
studies and simulations. Since POINTS is a "global 
astrometric instrument," all well-observed objects can and 
will be made to contribute to the stability of the reference 
frame used for each observation. In turn, each object that 
contributes to the reference-frame stability can and will 
be studied astrometrically and its motions modelled. A 
small number of observations of quasars will connect the 
POINTS reference frame to the best available candidate 
for an inertial frame. Table 1 lists the number of BL-Lac 
objects, quasars, and Seyfert galaxies of magnitude 15 or 
brighter from the catalog of Veron-Cetty and Veron. 11 
Finally, because of the - 90 deg separation of the target 
pair, for any target there are always several bright 
reference stars among the set of well studied objects. 


Our performance analyses assume the target to 
be an isolated, unresolved point source, i.e., one with no 
significant structure. In practice, this condition will not 
always be met. Star spots, flares, etc. will set the limit to 
the accuracy of the astrometric measurement of the center 
of mass of some targets. We may gain useful 
information about some stars from the motions of their 
centers of light. At 10 pc, the Sun has a magnitude of 5. 

A shift of its center of light by 0.1% of the Sun’s 
diameter, which is the largest one expects from sunspots, would cause an apparent astrometric shift of nearly 0.5 pas. 
For some active stars, the effect is much larger. The effect of a companion star is addressed in Section 4.2. 


Table 1. Number of reference objects of magnitude m or 
brighter 

m 

NUMBER OF OBJECTS 

BL-Lac Quasar Seyfert 

10 

1 

0 

0 

10.5 

1 

0 

0 

11 

6 

0 

0 

11.5 

9 

0 

0 

12 

23 

0 

0 

12.5 

38 

1 

1 

13 

61 

2 

1 

13.5 

99 

5 

3 

14 

166 

13 

5 

14.5 

266 

23 

7 

15 

381 

49 

15 


When an observation set has sufficient redundancy, it can be analyzed to yield a rigid frame; it serves to 
determine the angular separation of all pairs of observed stars, even those that were not observed simultaneously. The 
redundancy is measured by M, the ratio of the number of observations to the number of stars observed, which must be 
greater than about 3.5 to yield a rigid frame. With moderate redundancy (M=4.2). the uncertainty in the separation of all 
star pairs (a majority of which could not be observed directly because they are not separated by <= 90° ) is about eaual 
(on average) to the instrument measurement uncertainty. The grid is free of regional biases and may be further 
strengthened by additional data obtained when the grid stars are used as reference stars for additional science targets. 


If we use 300 grid stars plus a lew quasars and take M = 5, then the observation series requires about five days 
at the nominal rate of 350 observations per day. Such a series could be repeated four times a year to provide not only 
coordinates but proper motions and parallaxes for the stars. Our Monte Carlo covariance studies show that after ten years 
the coordinate uncertainties^ are -0.6 pas, the proper motion uncertainties are -0.2 pas per year, and the parallax 
uncertainties are -0.4 pas. 1213 Note that this parallax determination is a factor of 2 better than one would naively 
calculate from the coordinate uncertainties in a single series. The reason for this enhancement is that the 90 deg nominal 
angle of POINTS results in direct observations of "absolute parallax." 


Recall that in Section 2.4 w» noted that instrument biases can be determined from the analysis of an astrometric 
data set by virtue of the closure information it contains. The analysis would include a time-dependent bias model for 
each observing sequence such that the residual bias would be small compared to the measurement uncertainty. The 
appropriate form of the bias model will depend on a detailed understanding of bias mechanisms that we expect to acquire 
in later stages of design. Exploratory studies during the mission will be aimed at uncovering any bias mechanisms not 
identified before launch. The operational bias model will be augmented accordingly. All the recorded data from the 
mission may be re-analyzed at any time using a new bias model. 
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In other Monte Carlo covariance studies, we investigated the ten-year observing sequence with fewer 
observations. We found that observations can be deleted from the series by a variety of random or systematic procedures 
yielding an increase in the mean parameter-estimate uncertainty which depends principally on the square root of the total 
number of observations. Further, additional stars can be added to the observation sequence with a minimal number of 
observations (perhaps 20) per star. 

4 J2 Faint companions and immunity from biasing signals 

Some target stars will be undetected binaries. (The statistics of star densities on the sky indicate that accidental 
companions close enough and bright enough to have a biasing effect are too rare to be of interest except in particularly 
dense star fields, e.g., globular clusters.) We have shown that: (a) If the separation between the stars is small compared 
to the resolution of the interferometer baseline (VL - 0.05 arcsec), then the instrument tracks the center of light; (b) A 
well separated companion is ignored. For an unexpected companion one magnitude fainter than the target, at a 1 arcsec 
separation (in the interferometer plane) the biasing effect is small compared to the nominal measurement accuracy, and; 

(c) Between these extremes, the offset of the apparent position of the target from its true position varies as a sinusoid (of 
generally decreasing amplitude) in the star-pair separation. 

The above analysis assumed the star and its companion were modelled as a single source, and that they had the 
same temperature - a worst case. If the separation is comparable to or greater than the resolution of the interferometer 
baseline, or if the two stars have known different spectra, then their separation can be determined from the distortion of 
the channelled spectrum. With measurements of a known binary against several reference stars, the accuracy of the 
separation measurement is only slightly worse than would be obtained for the position of the dimmer star if the 
companion were absent. Equally important, the distortion of the channelled spectrum will reveal that the source is a 
previously unknown binary so that biased measurements can usually be avoided. Finally, we note that if the channelled 
spectra are preserved and if a previously undetected binary is revealed, the old data can be reanalyzed to removed the 
bias due to the companion. 

43 Data analysis, iterative and global nature 

Although the mission will support a large and diverse set of science objectives, we picture much of the data 
analysis as being done centrally because one worker's target is everybody's reference star. The effort required will 
depend on the degree to which we need to use spacecraft engineering data to help understand systematic errors. In any 
case, it is not a complex effort nor docs it require special computing facilities. If the data were available today, the 
central analysis could be done on a desk-top workstation. The effort would require a single analyst and would result in a 
reduced dataset that would be sent to the participating scientists for further analysis leading to publishable results. 

Early in the mission, the analysis will be limited to positions determined during observing periods of a few to a 
few tens of days. Of course, the best a priori proper motions and parallaxes will be included in the analysis model 
although these will have little effect. Timely analysis early in the mission will permit instrument problems to be detected 
without undue loss of data. At about one year, it will become possible to estimate the full set of five parameters 
(position [2], proper motion [2], and parallax [1]) for each star. The stability of the solution should increase considerably 

during the second year. 

After about two years of observation, the postfit residuals from the data reduction will be used to investigate 
irregularities in the motions of individual stars. Their apparent motions will be analyzed both for the intrinsic scientific 
interest and to improve the stability of the resulting reference frame. Initially, the position of each star will be 
determined with respect to the mean of the positions of all the stars. After the initial modeling of the irregularities in 
stellar proper motions, the position oi each star will be determined with respect to the nominal modelled reference 
positions of all the stars. Based on recent Monte Carlo covariance studies, we know that this iterative procedure 
converges quickly. The global motion of the frame will be measured against the quasars included in the set of observed 

objects. 
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4.4 Planetary system detection 

. If tU ? K are Jupiter-size bodies orbiting nearby stars, the astrometric wobble seen from the solar system will be 
orders of magnitude larger than the POINTS single-measurement uncertainty. However, the wobble will be 
istmguishable from proper motion only when the observations span a significant fraction of an orbital period We 

^f°™ Cd M !Th C M L miSSIOn SimUla “ 0nS With Single pIanetS around stars of the reference grid. Plants orbiting stars 
the grid should be an easier analysis problem because omitting them from the model will not distort the reference 

grid solution Preliminary studies confirm that multiple-planet systems will make greater demands on the observing 
system and the data analysis. (However, note that observing the solar system from outside, the main signature woufd be 
from Jupiter. Saturn s longer period produces a signal similar to proper motion, and the signals from Mars and Earth are 
“udy ‘ WUh Un,CSS ^ $yS,em " ^ VCry " eart > y - heavily-studied) {££££££?* 

, h In the ^ simulations, it became clear that knowing the right answer made the analyst’s task much easier To 
make the simulations more realistic, we developed a "double blind" methodology. 100 neaiby stars were selected from 

mffalwr? *? 3S h . osts for ,he ^bous planets for all simulations. A person other than the analyst would prepare 
™7*, file for the s,mulat, ° n program. The file contained a "seed" for the random number generator and ranges^r 
ass, radius eccentricity, and number of planets. Until a simulation was completed, no one knew how many planets 

simulation ** *** ‘ ey ° rbl,ed ’ and ,he analyst didn ’’ know what ranges of the parameters were used for the 

H SearChing for planets | s an iterative procedure. The analysis procedure identifies candidates for planet modelling 
a d op t‘ n J lzes Parameters to minimize the RMS error. Since this is not a linear system, these two faslrc interact initially 
srarTr^e ,ndud ? s on,y star coordinates, parallaxes, and proper motions. Postfit residuals for observations involving each 

a co t { freque r s - The ,argest signais flag ** »*« 

-.hp"!*i ® ut ’ since f star may he measured against one or more stars which have unmodelled companions power 

, S n0t a ' WayS COrTeSP ° nd '° 3 reaI Planet ' For «onately, the procedure is rebust; spuriouTZet ratals 

converge to near zero mass or otherwise reveal themselves as the iterations proceed. 

i on S,mUla r S h3d 3 , l00 ' star grid and 20 ' 50 P ,anels wi,h signatures ranging from 0.002 to 1 "Jove" c and periods 
of 1-20 years^ A mission lasted 10 years, and all observable star pairs were observed quarterly with the articulation range 

“ d wh^MW tvf raT* de,eC ' Cd W '' h greater ,ha " 90% P robabi,i 'y for signatures between 0.01 and 0.05 Jove 8 

defem liTgn^of Se " S ‘! ,V,ty dr0pped rapidly be,ow 0 01 Jove - A » ** orbital elements could be ’ 

d ^™ Cd U1 80% „ of ‘ he cases Wllh signatures larger than 0.01 Jove. Where not all orbital elements could be 
determmed, generally the observing geometry was edge-on or the orbital period was substantially longer than the mission 
length. Detections were reliable: there were no false alarms in the studies 8 

4*5 In-orbit testing and performance verification 

mr* three approaches 10 the in -° rbiI evaluation of the instrument. On the shortest time scale a few pairs of 

ihTmT ** C , Sen ^ S cal,brators - The y Wll > ^ observed regularly to provide a measure of short-term reproducibility of 
the measurements. Throughout the mission, the instrument bias will be estimated. The history of the estimated bias will 
be compared to pre-flight engineering estimates of its characteristics. Additional least-square^solutions for astrometric 
parameters and biases will be obtained with an enlarged bias set as a means of uncoveSg 

^rT,‘r, F n ay ’ le P ° Stfil rCSldualS fr ° m the analysis wi " 3 measure of performance ofUie instrument 
Astrometncally uninteresting, isolated stars will be important in understanding the instrument performance. 
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6 . APPENDIX: HIGHER ORDER RELATIONSHIP OF MEASUREMENTS TO STAR-PAIR SEPARATION. 


Each of the two POINTS starlight interferometers measures the offset 8 of a target star from the interferometer s 
optical axis. In the simplest description, the separation 6 of the targets is related to the separation <p of the optical axes of 
the starlight interferometers by those offsets 8 r where i = 1 , 2 . 


0 = <p - 8j + 8-, 


(Al) 


By design, 9 = n/2 + A, where the absolute value of the articulation angle A is up to 3 deg. Here we develop the next 
level of approximation to 6 , taking into consideration the effect of transverse instrument pointing errors and the 
differences between the orientations of the baselines and the corresponding pseudobaselines. 


Each POINTS starlight interferometer has a pair of fiducial blocks, which serve to connect the end points of a 
set of laser gauges. In particular, each laser gauge ends in a hollow retroreflecior on the fiducial block. The fiducia 
block design has the apices of the retroreflectors coinciding (to within a few microns) at the fiducial point. A fiduc 
point serves as a surrogate for the corresponding aperture; the pair of fiducial points m a single starlight interferometer 
define a pseudobaseline that by construction is very nearly the same as the real baseline, except for a translation, w 
mainly in the direction of the interferometer s optical axis. 


Let Pj (i = 1,2) be the pseudobaselines and Z Q the direction perpendicular to the "instrument plane” defined by 


Z 0 = P, x P 2 /|P,xP 2 | 


(A2) 


The corresponding interferometer axes (pseudo-optical axes) are 


A- = P: x Z n 


(A3) 


for i = 1, 2. The articulation angle <p is the angle between the unit vectors A l and A,, both of which are in the 
instrument plane. 

The real baseline B can be rotated from the pseudobaseline P as shown in Fig. Al, where C is the baseline 
rotation around the Z G axis and $ is the baseline skew toward the Z Q axis. For each starlight interferometer, we define an 
"interferometer plane." The interferometer axis Q is perpendicular to the true baseline. Since this leaves an 
unconstrained degree of freedom, we require the interferometer axis to be in the instrument plane. Then the 
interferometer plane is defined to contain both the interferometer axis and true baseline. The perpendicular to 

interferometer plane i has a direction Z t = Q, x B|. 

Shown in Fio Al is the target-axis offset. The target is at T, and is related to the coordinate frames by the 
following arc lengths' (1 - UD, 6 = QD, £ = DT, and v = VT. The otto is described as an in-plane displacement 6 . 
which is the quantity that the starlight interferometer measures, and e, the displacement of the target T from th 
interferometer plane, which must be measured by oiher means. (The subscript are dropped where u will cause no 
confusion.) The displacement of T from A can more simply be described by the spherical coordinates p = AV and v - 
VT. Below, p and v are found as functions of the small angles C and e * In turn ’ 6 1S found as a functl0n of ^ V' 
and v. 
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Figure Al. Spherical triangles. The target is at T, and is related to the coordinate frames by the following arc lengths: 
P = UD, S = QD, e = DT, and v = VT. 


In Triangle QUD, the required relations are: 


sin P = sin5 sin£ 

(A4) 

sina = tanp / tan!; 

(A5) 

cosk = sin!; cosa 

(A6) 
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where k = QDU and £ = DQU. In triangle ZqTD, the required relations are: 


sinv = cose sinP + sine cosp sinK 


(A7) 


siny * sine cosk / cosv 


(A8) 


We next replace the above expressions with the corresponding series expansions in the small quantities 6, and 
e, to obtain expressions for p and v. From Equ. A4: 


sinP = ; 


cosp = 



(A9) 


From Eqs, A5 and A9: 


cosa = 1 


2 


a 



From Eqs. A6 and A 10: 


COSK 




sin k = 1 


T 


From Eqs. A7, A9, and All: 


cosv * 1 


_ - e8£ ; 
2 


= e -if 


(A10) 


(All) 


(A12) 


From Eqs. A8. All, and A12: 


7 * 


We can now find p = AV = £ + a - y from Eqs. A10 and A13: 


(A 13) 


p = c + 8 - 



(A14) 


We next consider the angular separation 6 of the target stars and find an expression that is more accurate than 
Equ. Al. Figure A2 shows the geometry. We consider spherical triangle Zq T, T 2 , for which 

cosG = sinv, sinv, - cosv, cosv 2 sinu (A15) 


where u = q - rt/2. Note that u is about the same as A, the deviation from n/2 of the separation of the interferometer 
pointing directions. Thus, with the present instrument nominal specifications. |u| < 3 deg. We introduce p = 0 - T|, 
which we expect to be a very small quantity (i.e., measured in microarcseconds), and obtain 
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cos0 = - sin(u + p) = - sinu cosp - cosu sinp 
By combining Eqs. A 15 and A 16 and neglecting terms of order v 4 , we obtain 


P 


V 1 + v 2 


tan u 


cosu 


(A16) 


(A17) 


We are now in a position to find 0 - <p - +p 2 + pi° second order in small quantities. 


where 


0-(p = D+ E+ F 

(A 18) 

D = 8, - 8, 

(A19) 

E - C 2 - C, 

(A20) 

”> *> 

„ e,&,+(e7+e;)sinu 

, - - 

(A21) 
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In Equ. A18, the first term (D) is the principal interferomeiric measurement. The second term (E) is a simple, fixed bias 
related to the rotations of the baselines with respect to the pseudobaselines. This quantity is held fixed by the null servos 
of the FAM system. The third term (F) is a complex bias, further complicated by the realistic assumption that the 
measured value of e ( - will have a bias that we will need to estimate. It appears that by making a string of astrometric 
measurements with a range of values for E|, E 2 , and u. it should be possible to determine ijj, and the biases in E| and 
Ej. This is more readily seen if one assumes that several measurements are made of a single pair of stars, but the 
conclusion does not depend on that assumption. Finally, we note that the third order terms in 6 - 9 , which were excluded 
from Eqs. A19-21, are negligible. 
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ABSTRACT 

POINTS, an astrometric optical interferometer with a nominal measurement accuracy 
ot 5 microarcseconds for the angle between a pair of stars separated by about 90 deg, is 
presently under consideration by two divisions of NASA-OSSA. It will be a powerftil new 
multi-disciplinary tool for astronomical research. If chosen as the TOPS-1 (Toward Other 
Planetary Systems) instrument by the Solar-System Exploration Division, it will perform a 
definitive search for extra-solar planetary systems, either finding and characterizing a large 

° r Sh ° wing Aat ±Cy 3X6 far less numerous than now believed. If chosen as 
e AIM (Astrometric Interferometry Mission) by the Astrophysics Division, POINTS will 
open new areas of astrophysical research and change the nature of the questions being asked 
in some old areas. In either case, it will be the first of a new class of powerful instruments in 
space and will prove the technology for the larger members of that class to follow. Based on 

mdica ^ on of the observational needs of the two missions, we find that a single 
POINTS mission will meet the science objectives of both TOPS-1 and AIM. 

The instrument detects a dispersed fringe (channelled spectrum) and therefore can 
tolerate large pointing errors. In operation, the difficult problem of measuring the angular 
separation of widely spaced star pairs is reduced to two less difficult problems: that of 
measuring the angle between the two interferometers and that of measuring interferometrically 
the small offset of each star from the corresponding interferometer axis. The question of 
systematic error is the central theme of the instrument architecture and the data-analysis 
methods. Stable materials, precise thermal control, and continuous precise metrology are 
un amental to the design of the instrument. A preliminary version of the required picometer 
laser metrology has been demonstrated in the laboratory. Post-measurement detection and 
correction of time-dependent bias are the essential elements in data analysis. In that post- 
measurement analysis, individual star-pair separations are combined to determine both the 
relative positions of all observed stars and several instrument parameters including overall 
time-dependent measurement bias. The resulting stellar separation estimates are both global 
and bias-free at the level of the uncertainty in the reduced (i.e., combined and analyzed) 
measurements. 

1 INTRODUCTION 

Astronomical optical interferometry was first practiced by Michelson, who added a 
pair of translatable, fixed-angle siderostats to a large telescope so that he could measure the 
fringe visibility of stars as a function of interferometer baseline length. In 1920, he formally 
announced his determination of the diameter of Betelgeuse, to much public excitement. Each 
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of his measurements was a tour de force, and the technique fell into disuse for decades. (For 
a review of Michelson’s life and work, see Physics Today, 40(5), 1987.) More recently, new 
technologies have made possible construction of ground-based optical interferometers that 
could make routine observations of high scientific value. Early versions of these instruments 
are now working or near completion. Mariotti (1993) lists 17 ground-based interferometers, 
more than half of which are working. However, ground-based interferometers and telescopes 
jointly suffer two classes of restrictions: (1) Observations are made through a turbulent 
atmosphere with limited windows of visibility; and (2) The platform (Earth) rotates slowly, 
necessitating compensatory pointing, restricting the portion of the sky that can be observed at 
one time, and setting the time between observations of points well separated on the sky. 

These restrictions are removed when the observations are made from space. 

Eventually several optical interferometers can be expected in space. Support for this 
assertion comes from the report of the Astronomy and Astrophysics Survey Committee 
(AASC), established by the National Research Council under the chairmanship of John 
Bahcall to set the priorities for astronomy and astrophysics in the US during the 1990’s, and 
from the more detailed discussion of the future of optical interferometry contained in the 
report of the Interferometry Panel of the AASC (National Academy of Sciences 1991). 

Further support comes from a rapidly growing community of scientists and technologists 
interested in optical interferometry. We expect that most of the starlight interferometers in 
space will be imaging devices on a grand scale, with corresponding price tags. Because they 
represent a new and complex technology, these "greater observatory" class instruments are not 
likely to be flown without the benefit of a smaller, less expensive precursor instrument such 
as POINTS. 

POINTS, a space-based astrometric optical interferometer with a nominal measurement 
accuracy of 5 microarc seconds (pas), will be a powerful new multi-disciplinary tool for 
astronomical research (Reasenberg 1984; Reasenberg et al. 1988). The instrument includes 
two separate stellar interferometers that have their principal optical axes (nominal target 
directions) separated by <p, an adjustable angle of about 90 deg. (See Fig. 1.) Since the 
targets are within a few arcsec of their respective interferometer axes, off-axis distortions and 
the attendant biases are minimized. In each interferometer, the starlight leaving the 
beamsplitter is dispersed to form a "channelled spectrum," which gives the instrument an 
enhanced tolerance of pointing error, especially during target acquisition. A single 
measurement determines the angular separation of a pair of target stars. Each interferometer 
has a baseline 2 m long, and two afocal telescopes with 10:1 compression. The two 
interferometric subapertures of each starlight interferometer are 25 cm in diameter. Our 
studies of detectors and mirror coatings show that 20% of the photons entering the aperture 
will be detected. We find that the nominal 5 pas measurement uncertainty is reached after 
about one minute of observing a pair of mag 10 stars. The largest error component is the 
photon statistics of the starlight. Based on present models of slew, acquisition, and dam 
collecting times, we estimate that POINTS would observe about 350 star pairs per day if the 
targets were mag 10. 
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oln . , f*** the interfe rometers are mutually perpendicular, the instrument naturally performs 
global astrometry, which was first practiced by HIPPARCOS, and provides three advantages 
(1) Instrument bias can be routinely determined by 360 deg closure. TTius, for reasonable 
observing sequences, one obtains a global reference frame that is free of regional biases to the 
nominal accuracy of the data. (2) Parallax measurements are absolute, not relative. (3) For a 
gjven target star, the reference star is chosen from the great-circle band of sky 90 deg awav 
maximizing die probability of finding a suitable reference. If the relative rotation range of 

^r°, m u e ? r ° meterS 1S ±3 deg ’ ±1S observable ban d has an area of 2160 square degrees 
>5% of the sky) and can be expected to contain about 80 stars as bright as (visual) mag 5- 
1200 stars, mag 7.5; and 17,000 stars, mag 10. TTius, the observation time is not dominated 

ZtI P t 0tOn , Fate f dim reference stars - Each rcfercn ce star is a carefully studied target 
star, chosen from the reference grid described in Section 6.1. Thus the motions of the * 

reference stars are well modelled and do not corrupt the measurements of other stars. 

M „. , The ^?^y. rnen ‘. U ^ es two kinds of metrol °gy systems. The first is the Full Aperture 
Metrology (FAM), which measures the Optical Path Difference (OPD) in the starlight 
interferometers. In each interferometer, two laser gauges measure the starlight OPD and hold 
it fixed with respect to the "pseudobaseline" defined by a pair of "fiducial points " The 

^°" d Sy t St ® m . uses tbese flducial Points to measure the angle between the pseudobaselines of 
the two starlight interferometers. or 


f„ii n FAM ’ i y odulated laser J ight is injected backward into the starlight optical path and 
fully illuminates the optical elements that transfer the starlight. The first element that 
encounters the starlight has a low-amplitude phase-contrast holographic surface. This 



Figure 1. Proposed POINTS configuration. Left, stowed for launch. Right, operating 
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holographic optical element (HOE) diffracts about 2% of the FAM light, and the samples of 
FAM light from the two sides of the interferometer are brought together to produce an error 
signal that drives a null servo. In the process, as discussed below, the "fiducial blocks" are 
surveyed into position. Each fiducial block contains four incomplete retroreflectors, 
constructed such that their apices coincide with the fiducial point to within a few microns. 
Knowledge of the six distances defined by these four points determines the angle between the 
pseudobaselines of the two starlight interferometers. 

The second kind of metrology system consists of a set of six point-to-point laser 
gauges that measure those six distances. The angle <p between the real baselines is the 
quantity of interest; the angle between pseudobaselines is measured. Since these two angles 
differ by a small constant, the measurement of one yields a biased determination of the other. 
The bias is estimated from the closure information in the astrometric data, as discussed in 
Section 6. 

In section 2, there is a brief history of POINTS. Section 3 contains a description of 
the instrument. A central aspect of the instrument, the laser gauge is discussed in section 4. 
The design of the spacecraft is discussed briefly in Section 5. Mission operations, including 
the bias determination and correction, are discussed in Section 6. 

2 BRIEF HISTORY 

The POINTS concept can be traced to 1974 when I.I. Shapiro was asked to provide 
ideas for NASA missions in the distant future. Among the concepts he offered was optical 
interferometry in space, and his proposed applications included a second-order deflection test 
of general relativity. A few years later, after completing some preliminary analysis of a 
design that included two articulated starlight interferometers pointed in opposite directions, 
Reasenberg established a collaborative effort with K. Soosaar et al. at the C.S. Draper 
Laboratory, Cambridge, Massachusetts. 

For the Draper studies, each of two interferometers had a nominal baseline length of 
20 to 25 m and a pair of subapertures of from 0.5 to 1.0 m diameter. These studies yielded 
three major conclusions. (1) A right angle separation of the starlight interferometers would 
be an important improvement. This design to use widely separated targets predated the 
development of HIPPARCOS, which has become the first spacecraft to perform global 
astrometry. 1 (2) With modest thermal insulation, thermal time constants of major components 
would be hours to days. (3) Adding imaging capability would be exciting, but would add 
considerable complexity and cost and appeared inconsistent with the dimensional precision 
needed for the astrometry. The imaging capability required rotatable siderostat mirrors and 
the ability to observe with the baseline far from perpendicular to the line of sight, which 
required using an auxiliary delay line to compensate for the extra path to one of the mirrors. 


1 Kovalevsky (1980) defines five types of astrometry from very narrow field to global, depending on the scale of the 
target separation. He notes that HIPPARCOS is the first instrument able to do global astrometry at high accuracy. 
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In the late 70’s, the drawings of POINTS explicitly showed a long boom suDoortine 
die solarocculters that were needed for the light-deflection test of genend relativity ? * 

o'" t USC ° f P ° INTS f0r a Search for extra -s°lar Planets at the AAS/DPS 
Meeting #11 Qayton, Missouri, October 1979. (Reasenberg and Shapiro 1979) That year a 

spherical enclosure, as shown in Fig. 2, was added to reduce the variation of temperature wdth 
instrument orientation, and to decrease the radiation temperature anisotropy, whiJTwould 
distort the precision structure. (Reasenberg 1980, 1982) It was not realized until 1991 that 
the occulter intended for the deflection test could be made large enough to puT^e entire 
instrument in a permanent shadow, thus making the enclosure unnecessary. 

By 1 978, die baseline and aperture had shrunk to 10 m and 0.5 m, respectively 
(Reasenberg and Shapiro 1982) However, it was not until 1983 that they reacted thefr 
current nominal dimensions of 2 m and 0.25 m, respectively. This shrinLg reflected a 
Rowing realizanon of the close connection between the cost of a mission ate the size and a 
fornon the mass of the spacecraft. Of course, were there an opportunity to build the 
instrument with larger baselines or larger apertures, we would be pleased. However for now 
we seek to use technology m place of size to achieve the desired accuracy 


A{tm h At . thc be *“ l * of 1983 ’ Shapiro and Reasenberg moved to the Center for 

to S y por^S t < i rt y H hereafter i he P ° INTS lab0rat0ry was ^Wished. From 1985 
q ,. ’ ° NTS ” 5ceived support from the Smithsoman Institution through its Scholarly 

Studies Program. From 1985 to 1989, POINTS (along with a few other promts) ^eTvte 



Figure 2. POINTS with spherical enclosure, as envisioned in 1979. 
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modest funding through a CfA-wide grant from the Innovative Research Program (IRP) of 
NASA’s Office of Space Science and Applications (OSSA). For the three years starting early 
in 1989, POINTS was supported directly by the IRP. It was the work done under the 
Scholarly Studies Program and early IRP funding that developed the concept to a level that 
was attractive to the Solar System Exploration Division (SSED) of OSSA. The more 
substantial IRP support during the last three years permitted rapid progress in developing both 
the concept and critical technology. 

In 1988, the Solar-System Exploration Division (SSED) began to fund the 
development of POINTS under the Planetary Instrumentation Definition and Development 
Program (PIDDP). PIDDP has continued to fund the development as an instrument to search 
for and characterize extra-solar planetary systems, and perform additional solar-system 
research. That year, the SSED established the Planetary System Science Working Group 
(PSSWG) under the chairmanship of B. Burke (MIT). The PSSWG was to investigate 
opportunities for advancing planetary science by observing stellar environments and to 
recommend strategies to pursue these opportunities. An outgrowth of that SWG’s activities 
was the establishment by the SSED of the TOPS (Toward Other Planetary Systems) program. 
The program has three phases including TOPS-1, a space-based survey to find extra-solar 
planetary systems. The development of POINTS has been aimed at meeting the evolving 
requirements of the TOPS-1 mission, including that it meet the restrictions of a NASA 
Moderate Mission (for which the cost at launch + 30 days must be under $400M not 
including the cost of launch per se.) 

The TOPSSWG (formerly PSSWG) has recommended a set of requirements for a 
TOPS-1 search, which are included in the report "TOPS: Toward Other Planetary Systems" 
issued by the SSED during the last quarter of 1992. In discussing the TOPS-1 search for 
extra-solar planetary systems, the report states that (1) a negative result must be "a 
fundamental advance in our knowledge," and (2) "the search should encompass a sufficient 
population of target stars, ultimately exceeding 100 to 1000 stars of solar type, and ... extend 
to distances between 15 and 30 pc." (p. 18) Independently, the POINTS team has proposed a 
distribution of targets: 200 F, 400 G, 200 K, and 200 M stars not known to be in multiple-star 
systems plus a few hundred special stars, including binaries. Excluding those known to be 
binary and those for which there is insufficient information, we find that the faintest of the 
basic {F, G, K, M} set will be m b - 9.6. The most difficult astrometric signatures, based on 
a Uranus mass planet at Jovian distance (5.2 AU) from the star, range from 4pas for F stars 
(mag 6) to 50 pas for M stars (mag 9.6). In a POINTS mission concentrated on the search 
for extra-solar planetary systems, the targets would be bright and the required measurement 
correspondingly quick. A subset of easily observed targets (e.g., near) would be observed far 
more closely than the others to extract additional information in case they should have 
planets. 


In the past four years, there have been valuable contributions to POINTS by students 
working on the project. Two students have come to the project as fellows under the Summer 
Undergraduate Fellowships in Science program of the New England Consortium for Under- 
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graduate Science Education (NECUSE). That program, which places students in research 
situations outside their home institutions, has funded a major part of their stipends. Others 
a e worked as laboratory assistants and in undergraduate honors thesis work. 

at n Pa ?J Ck Morrisse y ^gan working with us through the NECUSE program while a junior 

POTm^ n H U ? 1VerSlty ' He 1 D Cr did hlS SCni0r h0n0rs thesis on some erTOr s °arces in the J 
NTS distance gauges. Patnck measured the polarization-dependent phase shift of the 

l3Ser g , aUgC endP ° im C ° mercubes ’ showin £ that the phase differences 
ere large, and that there were large variations from mirror to mirror (25° to 45°) but that 

they were very nearly constant with position on the mirror, and with rotation of the mirror 
about its normal Patnck did a significant amount of the work needed to commission the 
laser gauges, including careful identification of the causes of stray light in the laser 
modulator, and suppression of the stray light. He also showed that the sensitivity of the quad 

and T ng WaS Ve n ncarly COnStant 0Ver ** surfacc ’ 561 U P other optics, and built 

University^ lectron,cs - Patnck 1S now a graduate student in astronomy at Johns Hopkins 

Andrew Gerber worked with us for two summers, one as a NECUSE fellow. The first 
summer, he wrote a multi-channel data acquisition program for our laboratory PC His 

^ gram , aCqUirc u S data ’ displays h graphically, and allows changing parameters and writing 
cormnents mto the file without disturbing the data acquisition. The second summer AnZ w 
worked on an analytic ray tracing method, the forerunner of a system that we have used for 
the analysis of the POINTS high precision optics. Andrew wrore analy^aTrey ^cinf 
routines using algebra software running on machines of the KM PC family and identified the 
size of the expressions as the main problem that needed to be addressed. Llew !s now a 
fourth year student in physics at Yale. 1S now a 

on th, ST"? U " dergraduates Gordon p auth and Tom Killian have worked as lab assistants 
onto jncomciCT distance gauges. A large part of their work was building electronic equip- 
ment to our designs. They also worked with us testing the electronic systems that are ^ed the 
by the gauge, and analyzing the test results. 

3 INSTRUMENT DESCRIPTION AND CONCEPT 

3.1 Principles of Operation 

mnnnt T* 1 ! pnncipal elements of the POINTS instrument are two starlight interferometers, 
ounted at a nearly right angle, and a metrology system. The instrument determines 0, the 
angular separation between two widely separated stars, by measuring <p, the angle between the 
interferometers, and measuring independently 5, and the offsets of the target stars from 
their respective interferometer axes. With proper selection of target stars and a small 
adjustment capability m <p (nominally 87 < cp £ 93 deg), a pair of stars can be brought 
simultaneously to within less than one arcsec of their respective interferometer axes. Once a 
target star is m the field of an interferometer, 6 is measured through the analysis of the 
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dispersed fringe which forms a channelled spectrum. A central feature of the instrument is 
the real-time monitoring of the angle between the interferometers and the metrology along the 
starlight optical path, each of which uses a laser interferometer scheme based on currently or 
soon-to-be available technology. 

In each of the two interferometers, the afocal telescopes compress samples of the 
starlight, which are directed toward the fringe-forming and detecting assembly. At the exit 
ports of the beamsplitter, the light is dispersed and focused onto a pair of linear arrays of 
detectors. When the star is on axis (8 = 0), the signal at any given wavelength has equal 
intensity at the two beamsplitter exit ports. When the star is off axis (5*0), constructive 
interference for a given wavelength at one port is complemented by destructive interference at 
the other port At each port an alternating pattern of constructive and destructive interference 
is observed. The resulting complementary channelled spectra form the basis for determining 
5. Since the detectors in the array see narrow, contiguous portions of the optical spectrum, 
there are effectively a large number of narrow-band interferometers that collectively make use 
of all of the light. Because of their small bandwidths, these interferometers can function 
when the instrument is pointed several arcsec from the target. However, to keep the fringe 
visibility high, it is desirable to keep the pointing offset small compared to 8^ the Nyquist 
angle, at which there are two detector pixels per fringe. 

The channelled- spectrum approach has two distinct advantages over a system in which 
a single detector measures the white light fringe. In the latter system, the pointing error 
would have to be under 0.05 arcsec during observations (unless an active system of path 
compensation were included); white light fringe detection would require an oscillating mirror 
to modulate (say at sonic frequency) the OPD. The fringe visibility would be substantially 
reduced at 2 arcsec, making the initial acquisition of the fringe more difficult. Thus, the use 
of the dispersed fringe greatly simplifies the pointing system and removes the need for the 
motion of major parts (i.e., delay line) during an observation, which would add mass and 
complexity and reduce reliability. The second advantage of the dispersed fringe approach is 
that there is additional information available in the channelled spectrum. This information 
can be used to separate targets that are closely spaced and might otherwise be confused, such 
as members of a binary system. An instrument utilizing this technique exclusively was 
proposed by Massa and Endal (1987). 

There were two principal contenders for the coatings on the mirrors that transfer the 
starlight to the beamsplitter: silver and aluminum. Thus far, we have not considered complex 
coatings out of concern that they may have temperature-dependent phase shifts. Silver has 
higher reflectivity at the longer wavelengths, but is not useful at wavelengths shorter than 
about 0.4 micron. (The nominal bandpass for POINTS is from 0.9 microns to 0.25 microns.) 
In order to assess the relative merits of these two coatings and of various detectors, we 
developed a computer code that calculates the astrometric precision of one starlight 
interferometer by numerically integrating over a specified portion of the optical spectrum. 
This code can use various source spectra including a black-body spectrum (of specified 
temperature), a smoothed solar spectrum, and smoothed spectra of various types of stars. 
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fl< *! by ^ wavel ength-dependcnt reflectivity of the assumed mirrors 
e wavelength-dependent sensitivity of the assumed detectors to yield the detected 
spectrum. The study showed that to maximize the measurement rate for a TOPS-1 mission 
which would look at Sun-like stars, we should use silver coatings and neglect the ligTt 7 
wavelengths shorter than 0.4 micron; for an AIM mission we should use Aluminum Ld 

mirmrT ** W3Ve lengths ' However ’ we have more recently learned that silver 

^fnVrth Pr °, teCUVe COatingS ’ have a histor y of visible deterioration before launch 
therefore, the nominal mirror coating is aluminum. 

3.2 Detec tors and Instrument Pointing 

There are three questions concerning detectors and pointing: (1) How far from the 
target direction can the instrument point and still detect the target? (2) What is the effect on 
the mformauon rate of a pointing offset? (3) What is the liming nutgnitude of 2 
trument. As |5| is increased from zero, the number of fringes on the detector array 
increases, and the fringe visibility V decreases because each pixel averages over an 

V C S ° f t Mn fu ^ Ulf0 ™ ati0n Content of Ae data * monotonic with 

IniHw v. 7 ^ WhCre thCre 1S one fnnge P« P« cl - ^ information content falls 

rapidly with decreasing V near V = 1, so the visibility must be kept high. A prism, like a 

detector 00 ^ “r*® 0p P 0site sense )’ Provides a nonlinear dispersion so that, if the 

tSZXS? 7 d± ' ** nUm K ber ° f Pixds P 61 ^ of channelled spectrum 

* , ?" C S °! UUon 15 t0 combine a grating and a prism, as is being done for 

^ODtolifpo 1 !^ : merfer T ter 1990) - An alternative approach, which is 

adopted for POINTS is to use only a pnsm, allowing the pixel width to vary with the color 

thC ,ght l0SS ° f ** grating ' Below ’ we consider in tum the detectors, the 
pectrometer, limiting magnitude, and the fine pointing and isolation system. 

nh„t~i-!T ma c y thC nominal dete ctor had been a set of photon-counting avalanche 
photodiodes. Such detectors permit a fringe to be observed in the presence of a high fringe 

ate, /.e., when the instrument is rotating. However, as discussed below, we have found that 

a era frr ^ r, ff r, ntly Stable 0n 11,6 deteCtore * at “ detect™ such iL 

rrn U T Z U 1 C 0SS 0f visibmt y due to motion during the integration 
period. CCDs offer high quantum efficiency, and there is a wide base of experience of their 

wiVrrn* 6 ’ aS f wd ' as u co ; tinuin S development. Among the instruments that have flown 

n ° r t 7 f y “? devel °P ed 316 A e Wide Field/Planetary Camera I and 

H and the Space Telescope Imaging Spectrometer for Hubble Space Telescope, the Solid State 
Imaging camera aboard Galileo, and the imaging cameras and the Visual and Infrared 
apping pectrometer for Cassini. Since the spectrum is dispersed along the array, a tapered 
anti-reflection coating can be used, offering a further increase in quantum efficiency. With a 
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stabilized fringe integrated on the CCD detector, there is an automatic solution to the problem 
of data compression, which otherwise would require considerable on-board computation As 
presently envisioned, the CCD arrays will have 650 cells to cover an optical bandpass from 
0.25 to 0.9 pm. 

The spectrometers, one for each exit beam, must disperse the starlight to produce the 
channelled spectrum while maintaining high fringe visibility. Three primary effects govern 
resolution and thus fringe visibility: diffraction, aberration, and detector pixel size. In the 
dispersion direction, the diffraction-limited spot must be < 1/20 the spacing of the dispersed 

fringes in order to maintain high fringe visibility. Additionally, providing spatial resolution in 
the cross-dispersion direction allows simultaneous acquisition of data on multiple targets in 
crowded fields. With this capability, the scientific throughput on stars in globular clusters 
may be as much as two orders greater than on isolated targets. In the cross-dispersion 
direction, the limit to resolution is the diffraction limit of the subapertures. Both the 
requirement for high visibility and for resolution in the cross-dispersion direction dictate that 
the aberration-limited spot and the size of the pixels should be less than the Airy disk, i.e. the 
system should be diffraction-limited. 

The current design meets the above requirements. It employs a prism followed by an 
off-axis section of a paraboloidal mirror (f/40). The prism provides adequate dispersion with 
an apex angle as small as 10 deg. The nominal prism material is fused silica, because of the 
need for high transmission over the bandwidth from 0.9 pm to 0.25 pm. 

At one time, the limiting magnitude was set by the requirement that the interferometer 
collect enough photons to determine the fringe phase before the unmodelled part of the fringe 
drift smeared the fringe. The current design, however, calls for the fringes to be held 
sufficiently stable to allow on-chip integration, with duration limited only by the frequency 
with which instrument biases must be recalibrated. In this mode, limiting magnitude depends 
on dark current. We believe that present CCD’s impose a limit somewhere between 
magnitude 17 and 20. The two important factors for setting the limiting magnitude are 
detector temperature and total size. 

In the present design, the instrument is connected to the spacecraft scan platform by 
means of a soft, single-axis pivot with a free swing range of several tens of arcsecs. A 
flexural pivot that would serve well as the basis for the isolator has been made by the Perkin- 
Elmer division that is now Hughes-Danbury Optical Systems (HDOS). (The pivot was flown 
as part of the Apollo Telescope Mount on Skylab, 1973-1974. A. Wissinger, HDOS, private 


2 Consider an observation in which 10 2 * * * * 7 photons are collected in 100 batches by one interferometer. Without 

compression, two numbers (epoch and detector number, about ten digits total) must be transmitted to the ground for each 

photon. Without stabilization and on-chip integration, the fringe must be rotated and integrated in software, requiring several 

floating point operations per photon in addition to the real-time filter. With stabilization and on-chip integration, there is 

little computation required for fringe tracking, and the telecommunication load without explicit compression is 2 1(T numbers 

of five digits. It appears that a compression of at least 10 fold will be possible. 
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nrtim? UniCati0n - Eddy ! 979 ) 7716116 316 several suitable magnetic actuators to provide fine 
pointing correction to the isolated instrument. Tire instrument mounted on the pivot would 

have a torsional resonance between 1/20 and 1/5 Hz 3 . Given such a soft mounl! there would 
^considerable filtering at the detector sampling rate of 100 Hz, applicable to bright target 

The present scheme requires that at least one target be bright, say mag 10. Initially a 

.nTm^nT ™ USed 10 flnd the bri * hl *"*« ““I the approximate rotation rate of tire 
t^ ‘-t W ° U d ** used t0 reduce f* offset and rotation rate to below 

?! f f0F *“? m * s^atlight fringes. Once the fringes wen: found, they would serve as 

the reference for the fine pointing. On a mag 10 star, with a sample rate of 100 Hz the 
samples would have a stansncal uncertainty of -250 pas. However, the use of one of the 

!n h? a M m ' ttr5 ? ? e “ gUlar refercnce roontres real-time data analysis and the corresponding 
on-board computational capabthty. The 100 Hz sample rate probably limits to -10 Hz the 
unit-gam frequency of the control loop for the orientation of the instrument. 

3.3 Internal Measu rements and Systematic F.rmr 

The control of systematic error is central to achieving the stated instrument 
performance and the mission science objectives. We address this problem at three levels- m 

core^u^oT 18 T d thermal C ° ntr0l; (2) real ’ time m^ogy; and (3) the detection and 

of systematic error in conjunction with the global data analysis. The last of these 

s addressed in Section 3.4. The best materials fail by orders of magnitude to provide the 

fOPD^H 11 ^ 611510 "^ i^ blhty f ° r CaCh intcrfc rometer’s Optical Path Difference 

(OPD) and other critical dimensions, in the absence of other means of control. Stable 

StrUCtUnd dementS ° f Ae instrument servc to limit the dimensional changes 
at the metrology system must determine, and to control errors that are second order in * 

1 ntS ,n a plaCeS ' we ™ •» roly on material stability (over 

..." !,™ ) ' . '"smiment is designed so that such metering elements are small, and can be 
thermally isolated and regulated. 

For a 2-m baseline, the nominal 5-pas uncertainty corresponds to a displacement of 
one end of the interferometer toward the source by 0.5A = 50 pm (picometer = 10 12 m) 

Since similar displacements of internal optical elements are also important, the instrument 
requires real-time metrology of the entire starlight optical path at the few picometer level 

kHUr Tf n0t P ° Sean overwhelmin g Problem for two reasons. First, the precision 
£ needed only for averaging times between 1 and 100 minutes. Second, a slowly changing 
bias m the measurement is acceptable, as discussed in Section 3.4. 

. 1116 lnstru roent relies on two kinds of laser-driven optical interferometers to determine 

changes in critical dimensions. New concepts and technology have been required for the on- 


Wnf W wl!lch need t0 ^ K aWeS ’ e,C : CTOSSin8 lhe piVOt ' md ,hese W0U,d add 10 ** stiffness by an as yet undetermined 
mount, which is a contributor to the stated possible range of the frequency. 
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board measurements to picometer accuracy. The laser gauges per se, a new application for a 
holographic optical element (HOE), and the laser-gauge endpoint assemblies were invented 
during the development of POINTS. The laser gauges can be grouped according to whether 
they measure distances internal to a starlight interferometer or distances between starlight 
interferometers. Below, we consider the former. In the next section, we consider the latter 
and some related questions of instrument geometry (the relation between the baselines and the 
pseudobaselines) and measurement bias. 


The high-precision star position measurement is made with respect to the optical axis 
of the interferometer by measuring the difference in the optical paths from the target to the 
beamsplitter via the two sides of the interferometer. In turn, the position of this axis is 
determined (defined) by the positions of the optical elements used to transfer the starlight. 

The metrology system must determine, to about 10 pm overall accuracy, the average change 
in the starlight optical path differences (OPD) induced by all motions and distortions of all 
optical elements. Our approach is to use Full-Aperture Metrology (FAM), a novel technique 
discussed below. The roots of the FAM system in POINTS arc in the workshops held in 
1980 and 1981 by D. Black of NASA-ARC, and can be traced via D. Staelin of MIT to 
concepts developed for large DoD optical systems. FAM provides three significant 
advantages over conventional approaches. 

(a) FAM removes complexity. The usual metrology systems use a large number of 
laser gauges to determine the locations of the elements individually. From these 
measurements, the optical path through the system is computed. FAM directly measures the 
optical path through the system. 

(b) FAM measures the correct quantity. Because the metrology signal fully 
illuminates the surface of each optical element that determines the starlight OPD at the 
beamsplitter, the metrology OPD represents accurately the average starlight OPD through the 
system. 

(c) FAM provides the basis for an operational definition of the direction of the 
interferometer baseline. It measures the OPD with respect to a pair of fiducial points that 
define the pseudobaseline, which is located in front of each interferometer and which is held 
at a small fixed angle (nominally zero) to the real interferometer baselines. These fiducial 
points are used to determine <p, the angle between the two interferometers optical axes. (The 
small bias resulting from the angle between the pseudobaseline and the real baseline is 
determined routinely as part of the data analysis.) 


Figure 3 illustrates the preferred version of the technique. The key element is the set 
of primary mirrors which have shallow (phase contrast) zone plates on their surfaces: alternate 
zones are depressed about 100A. In this holographic optical element (HOE), the zones are 
approximately in the form of Newton’s rings; each zone has about the same total area. We 
have shown that the required Zone-Plate Mirror (ZPM) can be made holographically such that 
the diffractive spherical aberration is zero. Since the system functions essentially on axis, 
there are therefore no Seidel (/.<?., third order) aberrations in the diffractive focus. 
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In ^ AM ’ m °dulated laser light is injected backward into the starlight optical path and 
fully illuminates the optical elements that transfer the starlight The hologram diffracts about 
% of the FAM light, and the samples of FAM light from the two sides of the interferometer 
are rought together at the metrology beamsplitter to produce an error signal that drives a null 
servo. In the process, the fiducial blocks are surveyed into position. Each fiducial block 
contains four incomplete hollow comercube retroreflectors, constructed such that their apices 
coincide at the fiducial point to within a few microns. Knowledge of the six distances 
defined by these four fiducial points determines the angle between the pseudobaselines of the 
two starlight interferometers. 

For the flight hardware, we would manufacture the holographic optical element (HOE) 
by a method that would provide the required spatial dependence of diffractive efficiency The 
needed pattern would be found numerically and used to make a mask such as is used in the 
manufacture of integrated circuits. The mask could be written on a master plate by an "e- 
beam machine" or a "laser lathe." Each mirror would be given a thin (ca. 100A) coating of 
metal and contact printed using photoresist and metal etching. Finally, the optic would be 
given a reflective coating. This technique has been used for many years at Hughes-Danbury 
Optical Systems, Inc. (HDOS) in support of DoD projects. 

It is important that the uniformity of the thickness of the metallic layer on the optic be 
of controlled to about 1%, as the diffractive efficiency is proportional to zone thickness in the 
intended regime. We propose the following three-step process: (1) Normal sputtering. (2) 
Evaluation with a densitometer. The coating of interest will be about one optical depth. (3) 
Correction with Plasma Assisted Chemical Etching. With this technique, it would be possible 
to achieve a controlled non-uniformity to compensate for the (nearly Gaussian) profile of the 
laser light used to illuminate the ZPM. Senior engineers at HDOS have reviewed this plan 
and opined that it will work. Other approaches are being considered. 

The FAM technique requires two principal servos. The laser beams diffracted from 
the HOE’s in the two arms of the interferometer are recombined at a metrology beam splitter 
this is a Mach-Zehnder interferometer configuration, which is sinusoidally sensitive to the 
length difference from beam splitter to beam splitter via the two paths. A servo regulates this 
path difference to a constant value by moving the starlight beamsplitter assembly. A second 
laser beam injected into the metrology beam splitter parallel to the FAM beams similarly 
measures the difference in distances from that beam splitter to each of two comer cube retro- 
reflectors situated in the fiducial blocks; a second servo regulates that path difference to a 
constant value by moving the auxiliary beamsplitter assembly. With both servos working, 
there is a small (under 1 mm) and constant difference between the distances from the starlight 
beamsplitter to the fiducial points, which are the apices of the retroreflectors inside the fidu- 
cial blocks. These servos can have small bandwidths because the time scale for distortion is 
long compared to a second, and because we require that vibrations of the OPD be small 
enough that they need not be tracked. Since it is desirable to limit the contamination of the 
starlight by the laser signal ( e.g via scattering from the surfaces of the optical elements) a 
minimum of laser light is used in the FAM servo. 
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Figure 3. Interferometer optical paths, (left) Starlight and auxiliary 
the latter, a laser signal passes through the spatial filter A and is divided by the menology 
beam spli tter B to form beams that enter the fiducial blocks C and 1C . Within each fiducial 
btodMhe beam is deflected by the 45“ miiror, retroreflected by the hollow corner cube, and 
returned toward the metrology beamsplitter by the 45“ mirror The r^g agnals «e 
combined at the metrology beamsplitter and fall on detectors D and • ( E ) 
interferometer. A laser signal passes through the spatial min E and injection teamsph 
and is divided by the starlight beamsplitter G. The separated beams are reflected from the 
ternary minors H and H ' and secondary mirrors 1 and 1 ' to fully illuminate the pnnwy 
mirrors K and K 1 . The zone plates on the primaries diffract the signal to focal points in the 
^mHecondanes and .he athermai lenses collimate die diffracted light. FmaHy die 
sfgSk are reflected from the 45" mirror, are recombined at the metrology beamsplitter, and 

fall on detectors L and L . 
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1 qq-i\ I n e ^nF ematiVe f nd aPPr ° aCh t0 FAM < Reasenb ^ 1981; Reasenberg « al. 
1993), the HOE is on a large folding flat (i.e., fixed siderostat) which is the first element 

i'nterferom^ ^ Star lght< ' Phe two compressors (telescopes) in each of the starlight 
interferometers are mounted along the baseline, "looking” in opposite directions. In this 

angement, the fiducial block is located toward the target star from the folding flat, and is 
| “ ed "? the telescope secondly. The primary disadvantages to thts schenTare 
( ) the asymmetry of the geometry at the HOE, and (2) the need for an extra large optical 
element for each telescope - a total of four for the instrument. For a large scale instrument, 

1 * C ** a W !if ht WOuId clea [ ly ** im Portant; for a small instrument like POINTS, the case is 
less clear. There are two advantages to this scheme. (1) The HOE is external to and 

bv t^nthL r telescope. Thus the behavior of the diffractive optics is not complicated 
by the problems of telescope alignment. (2) The telescopes are "buried" deep inside the 
instrument, which makes them less susceptible to heating by stray light. 

, T* 1 ® fld “ cial blocks (Fig. 4) and the metrology beamsplitter assemblies pose the only 
. P roblem !dentified. To minimize changes of the size or shape of either of 

these, which would cause a corresponding changes in the bias of the metrolo^ system, they 
wi made of a stable material and kept in a thermally stable environment. For the fiducial 
*'n t £rf tlall: ' «”>*««* Premium ULE®; the beamsplitter material is still under 
m y ' j 1 ff “ ve separatl0n of metrology links within the fiducial blocks would be under 

hmiT ™ « Kmp ' ra ' Ure would h ' ld “Me *» 10 2 deg C on a time scale of a few 

hours. The effects of temperature fluctuations on a longer tune scale will be routinely 

removal during the analysis on the ground of the data. The thermal gradient will be reduced 
y putting an aluminum cylinder with 0.5 cm walls around the fiducial block We plan to 

““ rcqU ,' r ^'1 e T l " tab i' lty With acdve <* ^ outer skin of die fiducial 

S? M ™re, t0 10 d '8 C (although achieving a factor of ten better does not seem 
difficult) and a layer of insulation between the outer skin and the aluminum cylinder The 

IreldTT COn ^' ned WUh the h ' at capacitJ ' of flducial “ock and aluminum cylinder, will 
~ ° f ^ M 0m day - Usi "« Premium ULE®, which has a linear 
coefficient of thermal expansion |a| < 5 x 10-9/deg C in a range from 5 to 35 deg C, the 

Er^enBofmro beS i° T More important would be the effects of thermal 

reot.fr ° f ZrJ°- , d K? C “ V ” O'" 5 cm separation of metrology beams where they 
reflect from the fiducial block mirrored surfaces. Such gradients would produce a 
measurement bias under 1 pm, again within the required limit 


the -w £ n ? nu K fac k ture of flducial blocks ' the critical concern is for the alignment of 
the retrostnps, which are sections of hollow comer-cube retroreflectors. Each fiducial block 
has four retrostnps, each of which must be assembled onto the main body from three small 
pieces of glass Each retro has a symmetry axis that must point in approximately the correct 
direction and all four apices must be colocated to within a few microns. 


We now envision a manufacturing jig that manipulates and holds the three small 

P1< ?® S ° f f laSS in position on the main bod y- During assembly, the jig would permit the unit 
to be rotated on precision bearings, perhaps air bearings. The alignment procedure calls for 
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the retrostrip under assembly to form one part of an interferometer; laser light would be used. 
Under rotation, the changing fringe pattern would indicate misalignments of the three small 
pieces of glass. Algorithms have been developed for converting the interferometer fringe 
pattern into corrections to be made to the positions of the three small pieces of glass. The 
interferometric measurement is not expected to limit the accuracy of the alignment of the 
retrostrips. The precision bearings, which should have sub-micron runout, and the devices 
used to position the small pieces of glass will probably set the limit to the alignment. 

3.4 Determination of (P 

The angle <p between the baselines of the two interferometers (i.e., the instrument 
articulation) is determined by the measurements of the six distances among four fiducial 
points in the system. 

cosq>= (d 2 2 +d 3 2 -d 1 2 -d 4 2 ) /2b m b f 

Here dj are the four distances between the fiducial points of one interferometer and the 
fiducial points of the other and b m , b f are the baseline distances for each interferometer, i.e.. 



Figure 4. A fiducial block. There are four "retro-strips," each a slice of a truncated comercube 
retroreflector, and each to serve as the end point for a laser gauge. The apices of the comercubes 
coincide at the fiducial point in this computer-generated engineering sketch. 
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each is the separation of the fiducial points of one interferometer. Each of the six distances is 
measured by a laser gauge. (The laser gauges « are discussed in Section 4, where we 
ote that the techniques we have developed easily extend to an absolute distance gauge ) We 

-^m n r 8 r tC< iw e pr0blems associated with the laser gauges reading the distance correctly 
except for a fixed bias, as would be the case if an incremental gauge were used and the zero 

CrUde,y * ! meanS ^ - encode/on the tSlSon ^or 
from the known positions of bnght stars observed by startrackers mounted on the 

interferometer optical benches.) Under these conditions, there would be an additional five 

SyS " m W0U ' d 10 d ' ,arai "' d the ashomemr^feach 
nrae the laser gauges were turned on. Nominally, this would be at the start of the mission 

Stud,c ? show that ““mating these parameters poses no problem even if 
it were decided to estimate them ab initio as often as every day. 

, r .T* 1 * 5 J netrol °gy system measures the angle between the pseudo-baselines which are 

*"“! °f fldUCial P ° intS in the two starli g ht interferometers. The ’two (nearly) 
perpendicular pseudo-baselines define the "interferometer plane." For simple analyses we 

assume that each starlight interferometer lies in the interferometer plane. However we have 
considered the effects of the following small rotations of each of the starlight interferometers- 

(3) around tile ^ n * Cuh * t0 ** interferometer plane; (2) around the psfudo-baseline; and' 
^, lgh interferometer ’ s nominal optical axis. The first of these types of 
rotations yields a bias proportional to the difference of the rotation angles. The other four 

SSt? “tr"* in *“ fero "*«> d0 «* contribute to d mJL^l T 

fiist order, but do comnbute m second order. Thus the latter two types of rotations need to 
be ahgned to somewhat better than an aresec. In a related study, we arc considering how the 

£££££ ^erom™ ld ““ ^ » f * <**» « Aments of 

~ ss s"sis.“ ssstssa. 

mto the measurement of the angle <p. This bias is likely to be rime depends atTe 
mcroarcseclevel. lt is essential that we be able to determine and correct for the instrument 
HIPP A Rrn^H mth °* tbc introductio n of additional hardware. Because POINTS, like 
wh,n thJ ’ g astrometr y* bias determination and correction naturally occur 
rnnrHinof ob ^ xva,1 ^ ns combm ed m a least-squares estimate of the individual stellar 
oordinates (including proper motion and parallax), instrument model parameters and the 
expected brases. In particular, our covariance studies have shown that even w^outtiie 
mtroduction of a special observing sequence, it is possible to estimate simultaneously the 
stellar coordinates and tens of instrument bias parameters per day. The stellar coordinate 
estimate uncertainties increase by a factor of 1+0.003 N b , where N b is the number of bias 

SnTl ? Th CStimateC ! PCr ^ and M = 5 ‘ ^ ^ndancy factor, M, is defined in 
ThU ,!’ meff0 0 F biases 211(1 related erTors can he allowed to change on a time 
scale of hours without significantly degrading the performance of the instrument. This result 
allows us to relax the required stability at time-scales longer than 100 minutes of both the 


17 



laser used in the metrology and the rest of the instrument. This subject is also discussed in 
Section 6. 

4.0 Laser Gauges 

There are several standard means of using laser interferometry to measure distance. 

The European and American gravity-wave interferometer groups (Man et al. 1990, Shoemaker 
et al. 1991) have demonstrated Michelson laser gauges with more than 3 orders of magnitude 
higher precision than we need, but only for measurements of ultrasonic-frequency vibration. 
Salomon, Hils, and Hall (1988) have demonstrated DC performance about 5 orders of 
magnitude better than our requirement, but they used two-mirror cavity interferometers (not 
comer cubes) with very high quality mirrors. Commercial heterodyne distance gauges using 
corner cube endpoints (Hewlett Packard Model 5527 A, Zygo AXIOM 2/20 ) quote 

resolution no better than 1.25 nanometers, with systematic errors known to be much larger. 
Further, with the exception of a few complex and bulky devices, all of the laser gauges 
provide incremental distance from an initial epoch, not absolute distance. 

We have built a research prototype of one of the laser gauges we need using relatively 
crude interferometer optics and simple, inexpensive equipment. This gauge has met our 
no minal accuracy requirement for the flight instrument: 2 pm for averaging times from 1 to 
100 minutes. The "deviation" plotted in Fig. 5 is a root two-point variance directly analogous 
to the Allan variance which is used to characterize clock stability 4 . The laser gauge is the 
subject of a patent application filed 8 May 1992. 

Some of the POINTS laser gauges are intended to operate between comer-cube 
retroreflectors. If the laser beam makes an angle a with the line between the apices of the 
retros of such a gauge, the measurement is subject to a fractional error of l-cos(a). Although 
this error is proportional to a 2 , for few-picometer gauge accuracy we require half-arcsec 


4 One useful measure of gauge performance is an adaptation of the Allan variance used to characterize clock 
stability. For each averaging time x. the Allan variance is defined to be: 


m-l 


Ov ( t) = 


2 (m-l) 


E (y*-i ■ y*) 2 


k*l 


where y k are the frequency-offset averages scaled to the clock frequency. We use the term deviation to mean the square 
root of the analogous variance without scaling: 


m-l 


deviation = 


^ 2 (m-l) 


E (^.i - L *) ! 


k-i 


LOy (t) 


where L is the length being measured and L* are a set of such measurements averaged over interval X. This measures the 
predictability of drifting measurements as a function of the averaging and observation time scale X. 
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alignment stability. This can be achieved by a beam-position servo deriving its error signal 
from a quad cell behind a partially transmitting mirror in the retro. 

An enhancement to the performance of the laser gauge comes from operating the pair 
of ren-os as a cavity. TTus configuration permits the removal from the measured pafh of dl 
the glass (with its temperature-dependent contribution to optical path.) For the cavity to be 
optically stable, which is desirable but not essential for the gauge to work, there must be 
some focusing power in at least one of the mirrors that form the retroreflectors. In a recently 

nrn mP H C Hl S , y ’ 7? sho ^ ed that such cavities significantly reduce the effect of misalignment 
provided that the fringe-detecting photodiode has uniform sensitivity across its face. 

tn nffi.r W H ha I C conside l red the usc of two gauging systems for each measured distance, both 
dictanr#. ^ undanc y in the picometer-resolution incremental gauging and to enable absolute 
g ugmg. In particular, reducing the uncertainty in absolute distance to less than the 

ri:r cng r° uid comca,on «“* * ta—* 

frinee slinuin^TV recovtr ^ from a loss of and would obviate the question of 

fringe slipping. The repeatable and unambiguous resolution of the integer order of the 

interferometer requires much higher gauge performance than otherwise needed However 

even significantly less accuracy in the absolute distance measurement would be of value. ’ 

POINTS does not require extraordinary laser wavelength stability. Of the nominal 
measurement uncertainty, 25 picoradians, we assign 5 picoradians to the laser (in a 

?n^ a H U ?* SUm SCnSC) ; F ° rtUnately ’ *“* docs not correspondingly high laser stability 

In the determination of <p, we depend on the ratio of laser-determined distances and therefore 
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Figure 5. Deviation of the distance measured by the laboratory laser gauge. 
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laser frequency shifts cancel. (Secondary effects limit this cancellation at about the 10' 
level.) In the FAM system, the quantity being measured (with 10 pm accuracy) is the 
difference between two nominally equal paths, which should be under 1 mm, and laser 
frequency shifts would affect only this difference. Thus, the stability in either system need 
not be better than 10" 9 , and that stability is required only for time scales from 1 to 100 
minutes. 

Recently there has been a rapid evolution of the technology of lasers suitable for space 
missions of several years duration. There are at least two solid-state laser technologies which 
can meet the POINTS requirement: external-cavity diode lasers (Hamson and Mooradian 
1986), and diode-laser-pumped solid-state lasers (Byer 1988). The former have already 
demonstrated both drift over one hour and linewidth of under 100 kHz, which is a little better 
than required by POINTS. The latter have demonstrated sub-Hz linewidth when locked to a 
reference cavity. In either case, to achieve the required long term stability, it may be helpful 
to make the cavity of a stable material such as Premium ULE , which has a thermal 
expansion coefficient |a| 5 x 10' 9 / deg C; its temperature would be held stable over a 
period of hours to better than 10' 2 K. 

In diode-laser-pumped solid-state lasers, the pump diode can be made multiply 
redundant by coupling diodes to each of several facets of the lasing crystal. Even with 
multiple pump sources, diode lasers with a long life will be required. Units with 28 and 100 
mW output, amply sufficient for POINTS, have already been space qualified for two 
missions, with a specified Mean Time To Failure (in one of the missions) of 250,000 hrs (>25 
yrs). They are available from Spectra Diode Labs, Mountain View, CA. 

4.1 Spacecraft Gauge Design Considerations 

There were several architectural considerations in the design of laser gauges for the 
POINTS instrument. Both path-difference and point-to-point distance measurements are 
needed in the spacecraft. The fiducial blocks are designed around hollow corned cubes for 
three reasons. (1) The comer cube properly returns the beam even under rotations around its 
apex, as will occur during articulation. (2) A uniform thermal expansion of the glass does not 
change the measured distance. (3) The comer cube apex provides a well-defined endpoint of 
the measured length. In the design we use, the comer cubes form cavities, and there is no 
need to place glass in the optical path under measurement The light is injected and extracted 
through partially reflecting mirrors. 

Distance measurements need to be made to 10"^ A.. Such a measurement would be 
much more difficult if it needed to determine directly an arbitrary fraction of a wavelength 
than if it were locked to a convenient phase and used a null technique. We therefore 
extended our earlier successful work with the alternating frequency null gauge to the tracking 
frequency gauge described below. 
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It is advantageous to use modulation techniques to lock to a minimum or maximum in 
the transmitted intensity: DC drift can be quite small, and the modulation frequency may be 
chosen m a region of low photodiode and amplifier noise. For modulation, dithering the 
mterferometer length is a poor choice since the actuator that imposes it may add serious DC 
drift and synchronously-varying misalignment, and the optic that is shaken may be too 
delicate or heavy for rapid modulation. Thus the natural choice is frequency modulation 
( ever 1983). In turn, this demands that the interferometer output be sensitive to the laser 
frequency (implying an OPD of many centimeters in a Michelson). Thus the interferometer 
outputs must be periodic in both the length and the optical frequency. The period in optical 
frequency (c/OPD for a Michelson or c/2L for a cavity) is known as the free spectral range 
(FSR). To maintain the intensity minimum, the servo may either adjust the measured distance 
with an actuator or retune the laser to track the distance change. In the latter tracking- 

frequency approach, we measure the frequency offset needed to stay at the minimum; this is 
proportional to the distance change. 

4.2 La boratory Apparatus and Results 

Our prototype Michelson system incorporates both the distance actuator and tracking 
frequency approaches. Light from a single laser supplies two separate distance gauging 
systems, both systems measure the same OPD. The signal from the first gauge is used to 
servo the interferometer OPD to be a constant multiple of the laser wavelength X. The signal 
from the second gauge is used to regulate the frequency offset 5v required to reach an 
adjacent mterferometer null. Since the interferometer length is controlled to be an integer 
number of wavelengths, this frequency offset is expected to be constant. Any variations in 
that offset may be attributed to noise and systematic errors in one or both of the gauges and 
to the failure of the OPD servo (e.g., from insufficient gain) to remove the OPD fluctuations. 
This comparison technique, common in tests of clock and laser stability, allows us to test the 
gauge performance with picometer resolution despite micron-scale length variations due to air 
turbulence and the supporting structure. 

The first of the two gauge systems was based on an older design using 5 kHz toggling 
of the laser frequency. A portion of the light from a commercial Zeeman-stabilized HeNe 
laser is sent through a double-pass acousto-optic frequency shifter, which shifts the optical 
frequency by Bragg diffraction from a moving RF acoustic wave in a crystal. The RF drive 
for the acousto-optic modulator toggles at 5 kHz between two ciystal-controlled frequencies, 

32 and 40 MHz. The emerging laser beam likewise toggles at 5 kHz between two 
frequencies, 64 MHz and 80 MHz away from the frequency of the HeNe laser. The beam 
passes through a single-mode polarization-preserving optical fiber and a linear polarizer, and 
enters the Michelson interferometer. The interferometer has an OPD near 1 m, yielding an 
FSR of about 300 MHz. The physical structure is approximately temperature-compensated 
(-200 nm/deg C) such that the temperature coefficient of OPD due to air dominates. The two 
alternating frequencies of the injected laser beam land to either side of the intensity minimum 
to which we are locking. A photodiode detects the 5 kHz square-wave variations at the 
output of the Michelson interferometer due to detuning of the average optical frequency from 
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this minimum. A lock-in amplifier demodulates the signal at 5 kHz and feeds it to a servo 
filter. A high-voltage amplifier boosts this signal to drive a piezo-electric crystal, which 
pushes one of the interferometer endpoints. A second photodiode records the AM present on 
the beam at the input of the interferometer. This signal is demodulated in another lock-in 
amplifier and fed back to control the RF power that drives the acousto-optic device; this 
regulates the AM on the injected beam to essentially zero. The OPD measurements were 
made in air, and thus required considerable gain in the servo to remove die effects of air 
turbulence and temperature fluctuations, neither of which is a consideration in the 
performance of the gauge in the spacecraft. 

With the interferometer OPD thus fixed (as a multiple of the laser wavelength), we 
remeasure the OPD using the tracking frequency method. A second portion of the light from 
the HeNe laser passes through a second double-pass acousto-optic frequency shifter driven by 
a voltage-controlled oscillator (VCO) at about 138 MHz; this system shifts the laser frequency 
by a slowly varying amount. The beam passes through an electro-optic crystal, where it 
acquires phase modulation at 8 MHz with modulation index of order 1. This beam is 
combined with the first at a beamsplitter before entering the optical fiber, thus the two beams 
exit the fiber and enter the interferometer very precisely overlapped. The second beam (about 
300 MHz away from the first) interrogates the interferometer fringe one FSR away from the 
first The two forms of modulated light are detected by the same photodiodes, but the 8 MHz 
AM signals are separately amplified and sent to RF mixers for phase- sensitive detection. The 
demodulated signal from the second light modulator is filtered and sent to the frequency- 
control input of the VCO to keep the second beam tuned to die intensity minimum. The 
frequency of the VCO is monitored as an indication of the distance measurement errors, 
which may originate in either the alternating-frequency system or the tracking frequency 
system. The demodulated signal from the other photodiode is either monitored as an 
indication of the AM on the injected beam or used in a servo to control the 8 MHz AM 
(Wong and Hall 1985). 

Figure 5, which resulted from the recent work with this gauge, shows the deviation 
falling to around 2 pm for averaging times between 44 seconds and 6 hours. Note that this 
data record extended for more than 90 hours. This meets or exceeds the spacecraft 
requirements as presently understood. 

The apparatus presently contains only commercially-available equipment and some 
ordinary mechanical and electronic devices made at SAO. Systematic error sources that were 
controlled to the required level include electronic effects such as RF pickup and 
nonlinearities, and AM resulting from interferences within the optical fiber. However, 
because the beams from the two gauge systems overlapped so precisely in the interferometer, 
many systematic effects were suppressed; they include beam motion across the endpoint 
mirror surfaces, thermal gradients in the interferometer optics, differential cosine errors, and 
polarization variations. 
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5 SPACECRAFT DESIGN 


The POINTS spacecraft will contain the usual set of systems to support the 
instrument: telecommunications, attitude control, power, and computation. The spacecraft bus 
wiU be equipped with a single axis scan platform for the instrument, and may require an 
isolator to keep jitter from entering the instrument from the spacecraft The other two 
degrees of rotational freedom will be provided by rotating the bus in ineitial space. That 
rotation will be constrained by the requirement that the Sun shield keep the bus and 
instrument housing in shadow at all times. The bus’ degrees of rotational freedom are thus 
around the pivot that joins it to the Sun shield and (with the shield) around the direction to 
e un. A paper describing the spacecraft is in preparation by Schumaker et al. (1993). 


ithin the instrument, there is an articulation mechanism that sets <p. This mechanism 
must permit about 0.1 radian of motion and should have a resolution of order 0.1 arcsec: 2 

10 St l PS ‘ coursc ’ dunn 2 observations, the combination of the articulation pivot and drive 
must be stiff in all six degrees, especially in the rotation that changes cp. 


th c 7 ] ie solar cell array is mounted on a set of "flip-out" panels that form a skin around 

in C r”I!'S Cin ? end of thc s P ac ecraft See Fig. 1. As discussed below, this shield is important 
pernutting large coverage of the sky while not violating the requirement that the instrument 
enclosure be shielded at all times from the direct light of the Sun. 

The question of sky coverage is normally considered in the context of an instrument 

oArv?rc ing u e P ° inUng dircction and exclusion angles for the Sun, Earth, and Moon. For 
. , rs, the problem is more complicated because the instrument needs to look 

simultaneously in two directions about 90 deg apart. Thus it takes four angles to describe the 
phase space of the sky coverage. For the purpose of calculating sky coverage, we constrain 
the target separation to be 90 deg, reducing the number of independent angles to three. We 
may neglect both the Earth and Moon exclusions since we are willing to wait a few hours for 
any particular measurement. (The nominal orbit is circular at 100,000 km, which has a four 
cUy period.) Thus we are left considering only the spacecraft and the Sun. The line between 
them is an axis of symmetry. Rotation about that axis can be neglected in discussion and 
analysis of sky coverage; there are only two nontrivial angles. 


In our sky-coverage study, the first parameter was cos(y), where y is the angle between 
the anti-Sun direction and the first star, as seen from the spacecraft. The second parameter 
was co, the Sun-star l-star2 angle on the celestial sphere centered on the spacecraft, i.e., the 
angle between the planes defined by: (a) spacecraft. Sun, and starl; and (b) spacecraft,’ s tar 1, 
and star2. The study maps out the portion of phase space that is accessible. With this choice 
of parameters, the area in the "cos(y) - co plane" is proportional to the volume of phase space. 

• no ^ n ® mmal configuration, 83% of phase space is accessible at any one time, and stars 
m 99.1% of the sky can be observed, but some with a limited range of CD. (Stars in 98.5% of 
the sky can be observed with a cumulative range of co of at least 180 deg.) In all cases, the 
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stars or star pairs that cannot be seen at a particular time can be seen a few weeks earlier or 
later. 


The original work has been extended to address the relationship between the available 
volume of phase space and the size of the solar shield. Increasing the size of the shield 
increases the sky coverage for shields similar in size to the nominal, which has a 30 m 2 area. 
However, saturation of sky coverage is found for shields with diameters larger than about 7m. 

6 MISSION OPERATIONS 

Here we consider five aspects of mission operations: normal operation of the 
instrument, the effect of a faint companion, data analysis, the search for other planetary 
systems, and in-orbit testing. We do not discuss detailed target selection since there is an 
overabundance of interesting targets and the actual selection will depend strongly on which 
division of OSSA sponsors the mission. 

6.1 Normal Instrument Operations 

Our understanding of the characteristics of the instrument’s operation is based on a 
series of covariance studies and simulations. Since POINTS is a "global astrometric 
instrument," all well-observed objects can and will be made to contribute to the stability of 
the reference frame used for each observation. In turn, each object that contributes to the 
reference-frame stability can and will be studied astrometrically and its motions modelled. A 
small number of observations of quasars will connect the POINTS reference frame to the best 
available candidate for an inertial frame. Table 1 lists the number of BL-Lac objects, 
quasars, and Seyfert galaxies of magnitude 15 or brighter from the catalog of Veron-Cetty and 
Veron (1991). Finally, because of the •= 90 deg separation of the target pair, for any target, 
there is always a large number of bright reference stars among the set of well studied objects. 

Our performance analyses assume the target to be an isolated, unresolved point source, 
i.e., one with no significant structure. In practice, this condition will not always be met. Star 
spots, flares, etc. will set the limit to the accuracy of the astrometric measurement of the 
center of mass of some targets. We may gain useful information about some stars from the 
motions of their centers of light. At 10 pc, the Sun has a magnitude of 5. A shift of its 
center of light by 0.1% of the Sun’s diameter, which is the largest one expects from sunspots, 
would cause an apparent astrometric shift of nearly 0.5 pas. For some active stars, the effect 
is much larger. The effect of a companion star is addressed in Section 6.2. 

When an observation set has sufficient redundancy, it can be analyzed to yield a rigid 
frame; it serves to determine the angular separation of all pairs of observed stars, even those 
that were not observed simultaneously. The redundancy is measured by M, the ratio of the 
number of observations to the number of stars observed, which must be greater than about 
3.5 to yield a rigid frame. With moderate redundancy (M=4.2), the uncertainty in the 
separation of all star pairs (a majority of which could not be observed directly because they 
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are not separated by « 90°), is about 
equal (on average) to the instrument 
measurement uncertainty. The grid is 
free of regional biases and may be 
further strengthened by additional data 
obtained when the grid stars are used as 
reference stars for additional science 
targets. 

If we use 300 grid stars plus a 
few quasars and take M = 5, then the 
observation series requires about five 
days at the nominal rate of 350 
observations per day. Such a series 
could be repeated four times a year to 
provide not only coordinates but proper 
motions and parallaxes for the stars. 

Our Monte Carlo covariance studies 
show that after ten years the coordinate 
uncertainties are —0.6 pas, the proper 
motion uncertainties are -0.2 pas per 
year, and the parallax uncertainties are 
-0.4 pas (Reasenberg 1986; Reasenberg, 
1991.) Note that this parallax 
determination is a factor of 2 better than 


wuuiu uaivcxy calculate irom the 

coordinate uncertainties in a single series. The reason for this enhancement is that the 90 deg 
nominal angle of POINTS results in direct observations of "absolute parallax." 


Recall that in Section 3.4 we noted that instrument biases can be determined from the 
analysis of an astrometric data set by virtue of the closure information it contains. The 
analysis would include a time-dependent bias model for each observing sequence such that 
me residual bias would be small compared to the measurement uncertainty. The appropriate 
form of the bias model will depend on a detailed understanding of bias mechanisms that we 
expect to acquire in later stages of design. Exploratory studies during the mission will be 
aimed at uncovering bias mechanisms not identified before launch. The operational bias 
model will be augmented accordingly. 


In other Monte Carlo covariance studies, we investigated the ten-year observing 
sequence with fewer observations. We found that observations can be deleted from the series 
by a variety of random or systematic procedures yielding an increase in the mean 
parameter-estimate uncertainty which depends principally on the square root of the total 
number of observations. Further, additional stars can be added to the observation sequence 
with a minimal number of observations (perhaps 20) per star. 
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6.2 Faint companions and immunity from biasing signals 


Some target stars will be undetected binaries. (The statistics of star densities on the 
sky indicate that accidental companions close enough and bright enough to have a biasing 
effect are too rare to be of interest except in particularly dense star fields, e.g., globular 
clusters.) We have shown that: (a) If the separation between the stars is small compared to 
the resolution of the interferometer baseline (A/L — 0.05 arcsec), then the instrument tracks 
the center of light; (b) A well separated companion is ignored. For an unexpected companion 
one magnitude fainter than the target, at a 1 arcsec separation (in the interferometer plane) the 
biasing effect is small compared to the nominal measurement accuracy, and; (c) Between 
these extremes, the offset of the apparent position of the target from its true position varies as 
a sinusoid (of generally decreasing amplitude) in the star-pair separation. 

The above analysis assumed the star and its companion were modelled as a single 
source, and that they had the same temperature - a worst case. If the separation is 
comparable to or greater than the resolution of the interferometer baseline, or if the two stars 
have known different spectra, then their separation can be determined from the distortion of 
the channelled spectrum. With measurements of a known binary against several reference 
stars, the accuracy of the separation measurement is only slightly worse than would be 
obtained for the position of the dimmer star if the companion were absent. Equally 
important, the distortion of the channelled spectrum will reveal that the source is a previously 
unknown binary so that biased measurements can usually be avoided. Finally, we note that if 
the channelled spectra are preserved and if a previously undetected binary is revealed, the old 
data can be reanalyzed to removed the bias due to the companion. 

6.3 Data analysis, iterative and global nature 

Although the mission will support a large and diverse set of science objectives, we 
picture much of the data analysis as being done centrally because one worker’s target is 
everybody’s reference star. The effort required will depend on the degree to which we need 
to use spacecraft engineering data to help understand systematic errors. In any case, it is not 
a complex effort nor does it require special computing facilities. If the data were available 
today, the central analysis could be done on a desk-top workstation. The effort would require 
a single analyst and would result in a reduced dataset that would be sent to the participating 
scientists for further analysis leading to publishable results. 

Early in the mission, the analysis will be limited to positions determined during 
observing periods of a few to a few tens of days. Of course, the best o priori proper motions 
and parallaxes will be included in the analysis model although these will have little effect. 
Timely analysis early in the mission will permit instrument problems to be detected without 
undue loss of data. At about one year, it will become possible to estimate the full set of five 
parameters (position [2], proper motion [2], and parallax [1]) for each star. The stability of 
the solution should increase considerably during the second year. 
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After about two years of observation, the postfit residuals from the data reduction will 
be us«i to investigate irregularities in the motions of individual stars. Their apparent motions 
will be analyzed both for the intrinsic scientific interest and to improve the stability of the 
resulting reference frame. Initially, the position of each star will be determined with respect 

mean of '*■' P° su r s of a* *°rs. After the initial modeling of the tiregulariti« il 
stellar proper motions, the position of each star will be detennined with respect to the 
nominal modelled reference positions of all the stars. Based on recent Monte Carlo 
covanance studies, we know that this iterative procedure converges quickly The global 
motion of the frame will be controlled by the quasars included in the set of observed objects. 

6.4 Planetary system detection 


If there are Jupiter- size bodies orbiting nearby stars, the astrometric wobble seen from 
the solar system will be two orders of magnitude larger than the POINTS single-measurement 
uncertainty. However, the wobble will be distinguishable from proper motion only when the 
observations span a significant fraction of an orbital period. We performed Monte Carlo 
mission simulations with single planets around stars of the reference grid. Planets orbiting 
stare not in the grid should be an easier analysis problem because omitting them from the 
model will not distort the reference grid solution. Preliminary studies confirm that multiple- 
planet systems will make greater demands on the observing system and the data analysis 
This subject will require further study. y 


111 first simulations, it became clear that knowing the right answer made the 
a yst s task much easier. To make the simulations more realistic, we developed a "double 
blind methodology. A person other than the analyst would prepare an input file for the 
simulation program. The file contained a "seed" for the random number generator and ranges 
for mass, radius, eccentricity and number of planets. Until a simulation was completed, no 
one knew how many planets were included or which stars they orbited, and the analyst didn’t 
know what ranges of the parameters were used for the simulation. 

Searching for planets is an iterative procedure. The analysis procedure identifies 
candidates for planet modelling and optimizes parameters to minimize the RMS error. Since 
this is not a linear system, these two tasks interact Initially, the model includes only star 
coordinates, parallaxes, and proper motions. Postfit residuals for observations involving each 
star are examined for power at a comb of frequencies. The largest signals flag candidates for 
planet modelling in the next iteration. But, since a star may be measured against one or more 
stars which have unmodelled companions, power in the residuals does not always coirespond 
to a real planet. Fortunately, the procedure is robust; spurious planet models converge to near 
zero mass or otherwise reveal themselves as the iterations proceed. 


imulations had a 100-star grid and 20-50 planets with signatures ranging from 0.002 
to 1 Jove and periods of 1-20 years. A mission lasted 10 years, and all observable star pairs 
were observed quarterly with the articulation range set to yield M=5. Planets were detected 
with greater than 90% probability for signatures between 0.01 and 0.05 Jove, and with 100% 


27 



above 0.05 Jove; sensitivity dropped rapidly below 0.01 Jove. All six orbital elements could 
be determined in 80% of the cases with signatures larger than 0.01 Jove. Where not all 
orbital elements could be determined, generally edge-on observing geometry or orbital periods 
substantially longer than the mission length were found. Detections were reliable; there were 
no false alarms in the studies. 

6.5 In-orbit Testing and Performance Verification 

There are three approaches to the in-orbit evaluation of the instrument. On the 
shortest time scale, a few pairs of stars will be chosen as calibrators. They will be observed 
regularly to provide a measure of short-term reproducibility of the measurements. 

Throughout the mission, the instrument bias will be estimated. The histoiy of the estimated 
bias will be compared to pre-flight engineering estimates of its characteristics. Additional 
least-squares solutions for astrometric parameters and biases will be obtained with an enlarged 
bias set as a means of uncovering unexpected instrument behavior. Finally, the postfit 
residuals from the analysis will be a measure of the performance of the instrument. 
Astrometrically uninteresting, isolated stars will be important in understanding the instrument 
performance. 
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ABSTRACT 

Newcomb is a design concept for an astrometric optical interferometer with nominal single-measurement 
accuracy of 100 microseconds of arc (pas). In a three-year mission life, it will make scientifically interesting 
measurements of O-star, RR Lyrac. and Cepheid distances, establish a reference grid with internal consistency better than 
100 pas, and lay groundwork for the larger optical interferometers that are expected to produce a profusion of scientific 
results during the next century'. With an extended mission life, Newcomb could do a useful search for other planetary 
systems. 


The instrument is a highly simplified variant of POINTS. It has three (or four) interferometers stacked one 
above the other. All three (four) optical axes lie on a great circle, which is also the nominal direction of astrometric 
sensitivity. The second and third axes are separated from the first by fixed "observation angles" of 40.91 and 60.51 deg. 
The fourth axis would be at either 70.77 or 78.60 deg from the first. Each interferometer detects a dispersed fringe 
(channeled spectrum), which falls on a CCD detector array nominally 8k elements long and a small number of elements 
wide. With a nominal baseline length of 30 cm and optical passband from 0.9 to 0.3 microns, the Nyquist limit is 
reached by a star ±21 aremin from the optical axis. The instrument will be constructed of stable materials such as ULE 
glass, and have neither internal moving parts nor laser metrology. 

A reference frame can be constructed using stars located in a regular pattern on the sky. We start with the 60 
points that are vertices of the regular truncated icosahedron. The instrument axis separations are chosen as separation 
angles of the figure such that for each point there are four other points at the chosen angle. (There is no angle offering 
higher multiplicity.) To the original set of 60 points, we add two additional such sets by rotating the figure ±20.82 deg, 
which provides 40 interconnections among the three sets of 60 points. With such a grid of 180 stars, we have shown 
POINTS-like grid lock-up even with large Sun-exclusion angles (of up to 75 deg). In the covariance studies, we assumed 
nine quarterly observation series and five bias parameters per series per observation angle; we estimated position, proper 
motion, and parallax for each star. In an extension of the study, we eliminated stars at random from the set and found 
that there was stable behavior with as few as half the stars in the grid. Further, with 120 stars in the grid, the statistical 
uncertainty of star position increased by only 35% due to degeneracy (i.e., with constant number of observations.) 

As with POINTS, additional stars can be observed with respect to the grid. With three interferometers and a full 
grid and without exceeding the Nyquist limit, the region accessible for observing with respect to at least two grid stars is 
over 50% of the sky outside of the Sun exclusion zone. This observable region can be enlarged by either including more 
sets of 60 points in the reference star grid, or adding the fourth interferometer, or increasing the number of elements in 
the detector arrays, or by limiting the passband with the same number of elements. 
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1. INTRODUCTION AND HISTORY 

Newcomb is a design concept for an astrometric optical interferometer with nominal single-measurement 
accuracy of 100 microseconds of arc (pas). It is a low-cost derivative of POINTS (Reasenberg, et at. in these 
proceedings) intended for early demonstration of some interferometer technology and a reduced (from the POINTS goals) 
but significant scientific program in addition to the principal goal of establishing a high precision reference grid. 

The concept was first documented in a letter to M.S. Kaplan (NASA, Astrophysics) from R.D. Reasenberg on 11 
December 1992, describing a suggestion made by the latter in response to a concern raised by K.J. Johnston during a 
phone conversation of 10 December. There quickly ensued a collaboration between SAO and NRL to develop Newcomb 
for the Navy’s Space Test Program (STP). The plan was to use a standard spacecraft bus such as STEP (made by TRW), 
which could be put into Earth orbit by a launcher commonly used by the STP. The instrument would be kept simple and 
the mission, low cost and of short duration for an astrometric investigation, say 27 months. 

No detailed study has yet been made of systematic error in Newcomb. However, the plan to dispense with the 
metrology is plausible. If the instrument is constructed of solid ULE with expansion coefficient 1 x 10' 8 /K (in a worst 
case we may assume that this varies by 100% across the instrument) and thermal expansion of the 0.3 m dimensions must 
contribute no more than 30 microarcsec to the instrument bias, then the instrument average temperature and temperature 
gradients must be held constant to within -0.015 K over the recalibration timescale (several hours). Our thermal studies 
of POINTS have shown that this level of control is likely to be achievable, even accounting for the variation in thermal 
input from the warm spacecraft. 


2. INSTRUMENT DESIGN 

The nominal instrument comprises three or four Michelson stellar interferometers, each with a dispersed-fringe 
detection system. Each interferometer baseline will be parallel to the instrument’s "principal plane." As we picture the 
construction, there will be a stack of parallel plates with an interferometer assembled between each pair of adjacent 
plates. The nominal baseline length is 30 cm and the nominal aperture is 5 cm. Figure 1 shows the optical plan for one 
of the interferometers. 

Starlight from folding flats M, and M,' in each interferometer combines at beamsplitter BS. A beam from each 
beamsplitter exit port is directed into a spectrometer. The usual beam -compressing telescopes are absent, which reduces 
the number of components, the light loss in reflections, and beamwalk on Mj and M,'. The spectrometers employ prisms 
instead of gratings because of the light loss of the latter, which would be especially severe over the Newcomb optical 
passband (nominally 0.9 to 0.3 pm). The high loss would have contributions from the blazing, which could not be 
optimum over the whole passband. and from the filters needed to defeat the response at undesired grating orders. Given 
specifications for the prism, focussing optics, and detector, the range of optical frequencies 8v p corresponding to each 
pixel may be calculated. The blur due to the combined effects of diffraction and geometrical aberration expressed in 
terms of optical frequency must be less than Sv p at most or all pixels. This resolution can be provided by a prism with a 
68 deg apex angle. 

The detector has 8192 pixels of 7.5 pm width for a total length of 6 cm, and the prism deviation range is 4.1 
deg. The camera focal length of 85 cm matches the prism deviation range to the length of the detector. A single mirror 
design would provide the greatest throughput, but would require excessive length and would have unacceptable 
geometrical aberration over the range of angles produced by the prism. A two-mirror design was adopted. The camera is 
similar to the standard Ritchey-Chretien telescope, which is a Cassegrain with conic constants adjusted to cancel coma in 
leading order. 1 Dimensional stability of the camera will not be as much of a problem in this application as in ordinary 
telescopes as both the mirrors and their spacers can be made of ULE glass. A secondary magnification factor of 3 was 
chosen. A design with higher magnification would be more compact, but might have larger aberration and tighter 
tolerances. The conic constants were further optimized using the numerical ray trace program Zemax [FocuSoft, Inc., 
Pleasanton, CA], 
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As shown, the camera has a geometrical blur a factor of several greater than the diffractive blur at some 
wavelengths. Some of the options for further improvement are: additional aspheric terms in the shape of the mirrors (the 
present mirrors are hyperboloids of revolution about a common center line), an aspheric corrector plate, and reducing the 
deviation range (trading increased diffractive blur for reduced geometric). The optimization has not yet been performed, 
as it involves the dependence on optical frequency of all of the following quantities: the dispersion (i.e., the deviation 
angle per unit optical frequency), the diffractive blur, the photon rate from a typical star, the angular information in a 
single photon, and the photon detection probability. 

The starlight optics operate at non-normal incidence, and most mirrors exhibit a difference in phase delay for 
light polarized perpendicular and parallel to the plane of incidence (s and p polarizations). This implies that fringes in s 
and p polarized light will be formed in different places on the detector. The detector will record the sum of the two 
fringe patterns, so will see fringes of lowered visibility. The reflection phase difference will be minimized to the extent 
possible through coating design, but it might be necessary to have polarizing beamsplitters, and separate detectors for the 
two polarizations. 

Each of Newcomb’s 3 or 4 interferometers requires two of these spectrometers, one at each beamsplitter exit 
port. Figure 1 shows how two spectrometers might be arranged in relation to the interferometer optics. To make the 
design more compact, a secondary magnification factor higher than 3 could be used, or it might be advantageous to 
rearrange the components, reversing one or both prisms, and placing one or both secondaries on the side towards the star 
of the line from Mj to BS, and from Mj' to BS. 

3. THE STAR-GRID PROBLEM 


The mission concept for Newcomb relies on a reference grid of stars and a redundant set of intra-grid 
observations that permit the grid star positions to form a tightly connected system. This aspect of a Newcomb mission 
resembles a POINTS mission 2,3,4 for which the critical quantity is the redundancy factor: 


number of observations 
number of stars Jn the grid 


From our studies of POINTS, we know that when M > M 0 , where M 0 is about 3.5, the grid "locks up" such that the 
angle between any pair of stars can be determined by the analysis of the data, even when the pair of interest cannot be 
directly observed. For POINTS, we normally assume M = 5 in studies of the mission. In the nominal POINTS 
architecture, the angle between a pair of stars to be observed must be<t> = 4> 0 +A ’l A l- A o’ w here <t> 0 * 90 deg and 
A c is the one sided articulation range. If POINTS observes all of the star pairs within its articulation range, then 


A o sin <f> 0 

where N* is the number of stars in the grid. For POINTS, with the nominal parameters of Af=5, A 0 =3 deg, and <|>=90 
deg, we find that jV*=190. 

For Newcomb, if we were to have two interferometers at <|> o =90 deg, and if we assumed and A 0 = 21 min., 
then we would find that N* = 1629. If the instrument can make D observations per day, the time needed to make a 
complete set of grid-star observations is MxN*/D. Since D is of order 200, it would require 41 days to complete one 
observation cycle of the grid. This is clearly not acceptable. 

3.1. A sparse and regular star grid 

An alternative to the random grid of stars is a regular pattern that defines special angles, which are each the 
exact separations of a large number of star pairs. Such a "crystal on the sky" could be based on one of the semi-regular 
polyhedra. We selected one of the biggest of these, the truncated regular icosahedron (a.k.a. buckyball) as the starting 
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point. We have not yet tested several other semi-regular 
polyhedra that might form the basis of a reference grid. 

These include the snub icosidodecahedron, 
rhombicosidodecahedron, truncated dodecahedron, and 
truncated icosidodecahedron. Some of the dual polyhedra, 
which have vertices above the centers of the faces of the 
semi-regular polyhedra, also look promising. 

The truncated regular icosahedron has 60 vertices, 90 
sides, and 32 surfaces: 12 pentagons plus 20 hexagons. It is 
the pattern on the familiar soccer ball and the geometry of the 
nearly spherical molecule known as a buckyball. Each 
vertex can be mapped into each other vertex by a (possibly 
improper) rotation. Of the 59 angles from a given vertex to 
the other vertices, the following repeat four times: 40.91, 

60.51, 70.77, 78.60, 101.40, 109.23, 119.49, and 139.09 deg. 

No angle repeats more than four times. Observation angles 
near 90 deg provide maximum sky area in which to find 
target stars, but smaller observation angles are expected to be 
required for solar-glare isolation. A measurement grid must 
meet three increasingly stringent requirements: (1) it must 
lock up when a complete set of data are analyzed, (2) it must 
remain locked up when biases are estimated, and (3) it must 
be robust against deleting a moderate fraction of the stars. 

Grids were tested using the POINTS simulation program 
[TM81-6]. In a typical simulation, all observable star pairs 
were measured each quarter year for 9 quarters. Stars within a cone of adjustable size near the Sun were not observed. 
Five star parameters (two positions, parallax, and two proper motions) were estimated, as well as a variable number of 
Fourier bias parameters for each observation angle. Star positions were given a weak a priori estimate to break the 
overall rotational degeneracy, and all measurements were assumed to have 100 pas precision. 

A simple measure of grid robustness is the distribution of log 10 (inter-star angle uncertainty) for all star pairs. 

We prefer the mean of this to be somewhat smaller than log 10 (single measurement uncertainty), and that the distribution 
have few if any outliers. The 60-star grid is rigid with a single measurement angle of either 40.91 or 60.51 deg, but if 
even a single bias parameter is estimated each quarter, the bias uncertainty is large and the grid falls apart. Since M is 
only 2.0, the fragility of the 60-star grid is not surprising. 

A more robust system was assembled by using three sets of 60 stars and two measurement angles, 40.91 and 
60.51 deg. The second and third sets of stars are rotated from the first around an axis through the centers of opposing 
pentagonal faces by ± 20.82 deg to yield M = 4.22. (No intra-grid observation is possible at the 19.60 deg difference 
angle, but this angle can be used to observe non-grid stars.) The grid is robust against deletion of roughly one third of 
the stars, even if biases are estimated and near-Sun observations are excluded. Figure 2 shows the behavior of the 180- 
star grid, with 60 randomly-selected stars deleted, which reduces M to 2.85. Mean inter-star angle uncertainty increases 
slowly as the glare cone angle is increased to nearly 90° (keeping constant the total number of observations). Actually, 
insensitivity to solar glare was a general feature of all grids tested. If the grid were rigid, adding a glare exclusion cone 
merely increased slightly the uncertainty of the inter-star angle unless the cone angle was >90 deg, in which case some 
pairs were never observable and the grid unlocked. 

Sensitivity to missing stars was tested by making multiple runs with randomly selected stars deleted from the 
grid. When the grid is incomplete, some stars may become "orphans," that is. they do not have the minimum of three 
reference stars required for them to contribute to the information that locks the rest of the grid. Figure 3 shows scatter 
plots of uncertainties in inter-star angle for orphan, normal, and all pairs as a function of the number of stars in the grid. 
In selecting a set of real stars, we are free to try a large number of cases and pick the most favorable. Therefore, we take 



Figure 2. Degradation of star coordinate uncertainty with 
increasing solar glare exclusion angle. Dashed curve 
assumes that with fewer observations possible, 
correspondingly more time can be spent on each one. 
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Delete random points from 180 star grid 



Figure 3. Results from simulations in which randomly selected stars are deleted from 180-star Newcomb grid. 20 
runs were made for each grid size, (a) Distribution of numbers of orphan pairs encountered versus number of stars 
in grid. Numbers shown in the plot are the number of cases (max of 20) for which a particular number of orphaned 
stars was found, (b-d) Uncertainties in inter-star angles for orphan, normal, and all pairs. (An orphan pair is one 
which includes an orphan star. Orphan stars were not excluded when counting reference stars to determine whether 
another star is an orphan.) 
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the bottom of each vertical distribution as indicative of the 
results we expect from a formal search process. It appears 
that at least one third of the grid stars can be deleted without 
greatly increasing the angle uncertainties, more with selected 
matchings of the grid to the sky. 

3.2. Finding stars in the patches 

In the preceding section, we were concerned with the 
star grid as an abstract construct. Here we address the 
characteristics of a star grid based on real stars. When we 
select specific stars to form a grid on the sky, we require that 
both interferometers work within the Nyquist limit for all 
intra-grid observations. The simplest way of assuring this is 
to require that each grid star be within the Nyquist angle of 
the nominal position, i.e., to require that each selected star be 
within a circular patch with a radius equal to the Nyquist 
angle, and centered on the nominal position. Ideally, stars 
should be more accurately positioned than this to gel a higher 
information rate, but if we are unlucky, it may sometimes be 
necessary to observe well off axis with both interferometers. 
We can estimate how faint the grid stars might be from the 
average density of stars and the area of sky seen by each 
interferometer. If each star must be within 21 arcmin of the 
nominal position, then each grid star must be located within a 
circular patch of area 0.38 sq. deg. 

If the density of stars is p, and the sky area within 
the Nyquist angle is A, then the average number of stars in 
a patch is pA. Assuming Poisson statistics, the probability 
that a patch is empty is P 0 =e' pA . and the probability that a 
patch is occupied is 1 - Since the number of occupied 
patches has a Binomial distribution, the probability of 
having exactly m out of n patches occupied is 

C(n,m) (1 -e _pA ) m (e -P A ) n -™ (3) 

where C(n,m) is the number of combinations of n things 
taken m at a time. The probability of k or more occupied 
patches is obtained by summing Eq. 3 over the range 
m=A' to n. The result of this sum is shown in Fig. 4; we 
want to operate at P 0 <l/3, which implies pA>l.l. Table I 
gives the average density of stars on the sky as a function 
of visual magnitude. 5 Note that the star density is roughly 
a factor of two higher than tabulated near the galactic 
plane, and a factor of two lower near the galactic pole. 

From the table, we see that if stars were uniformly 
distributed, the grid stars would only have tc be as faint as 
m=9. The uneven distribution will require accepting fainter 
stars. However, as previously noted, we have the option to 



Figure 4. Probability that 120 or more of 180 possible 
patches of sky contain a grid star as a function of the 
probability that any single patch is empty. Insert is a 
blowup of the high-probability region on a linear scale. 


Table I - Density of stars 

m vis 

Stars/ 

0.38 sq. deg 

Expected no. of 
occupied patches 

5 

0.015 

2.7 

6 

0.046 

8.1 

7 

0.13 

22 

8 

0.38 

57 

9 

1.1 

120 

10 

3.1 

172 

11 

8.4 

180 

12 

22 

180 

13 

56 

180 

14 

138 

180 
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pick the best orientation of the set of patches on the sky and, by Monte Carlo methods, we can have a large number of 
cases from which to make the selection. 


33. Sky coverage 


A multitude of small patches of sky are observable relative to the reference grid. The centers of these patches 
are at P k which are the solutions of 


fV = cose ro 
Fy P k * cos 0 n 


(4) 


where are unit vectors pointing at appropriate combinations of grid stars and 0 m> 0 n are all possible combinations of 
observation angles. Physically, each measurement angle generates a band, bounded by a pair of small circles, 
surrounding each grid star, and the intersections of these with the bands of another grid star are observable patches. The 
comers of the approximately rhomboidal observable patches could be found by solving Eqs. 4 with 0, offset by the 
maximum off-axis angle. We took another approach. Sections of the sky were covered with a 0.1x0. 1° lattice and the 
number of references stars visible from the center of each lattice box was tabulated. (Because of the symmetry of the 
180-star grid, the study region needed only cover 36 deg of longitude and from the equator to one pole.) 


Figures 5a and 5b show the fraction of the sky visible against the nominal 180-star grid as a function of the 
inter-star angle tolerance. More than 30% of the sky can be observed with respect to at least the minimum of three 
references stars that are needed to estimate star-specific biases. Additional coverage could be obtained with additional 
60-star sub-grids. It might also be useful to put grid stars in the patches that can be seen from at least six of the grid 
stars. (Starting with a 60- or 120- star grid and adding stars at high-redundancy patches is an alternative way of building 
a grid, which has not yet been explored.) 


Visibility with 1 80 -star grid 




Inter-star angle tolerance (arcmln) 


(a) 


(b) 


Figure 5. Fraction of sky observable against the 180-star Newcomb grid. The inter-star angle is 19.6, 40.91, or 
60.51 deg, with a tolerance of ± 21 arcmin. Curves are labeled with N, (a) with exactly N reference stars and (b) with 
N or more reference stars. (Note that the 180 patches containing reference stars cover about 1/6% of the sky.) 
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4. APPLICATIONS 

Newcomb would establish a precise reference frame that would be accessible to optical and infrared sensors. 

This frame would be useful to the Navy for navigation. Hipparcos will yield a grid that is good for a short time, but 
proper motion uncertainty from this short-duration mission will quickly degrade that grid. Newcomb and Hipparcos 
would be complementary missions for determining proper motion. Repeated and redundant observations would insure 
that the program was robust and lower the position (and thus the parallax and proper motion) uncertainty well below the 
nominal single-observation uncertainty. 

Newcomb will also provide a direct link between high-precision optical astrometry and the present radio 
reference frame as discussed by Fey et a! 6 , and in the papers they cite. Observations of only a few radio quasars are 
needed to fix the relative rotation between the current radio reference frame and the optical reference frame which will be 
developed by Newcomb and Hipparcos. Newcomb will thus resolve current difficulties in relating high angular resolution 
observations at radio and optical wavelengths. 


Newcomb could see enough of the sky to do some interesting science. A preliminary study suggests that the 
instrument's limiting magnitude would be about 15. However, this analysis depends on instrument parameters that are 
not yet fixed and are expected to be set during a series of trade studies. The Space Interferometry Science Working 
Group (SISWG) developed a Strawman Science Program for the Astrometric Interferometry Mission 7 as a step toward 
evaluating the astrophysical capability of candidate missions. Of the goals presented in the "Strawman,’’ five could also 
be met by Newcomb: (1) 19 known Cepheids have parallaxes between 200 and 1000 pas and V magnitudes of less than 
10; a 5% distance measurement would be useful for refining the cosmic distance scale. (2) Absolute magnitudes 
of O stars are uncertain because none is close enough to allow a trigonometric parallax measurement from the ground. 

O stars are bright targets. V=4 to 6 at 1 to 2 kpc; 25 to 50 pas parallax measurements are needed. (3) Estimates of the 
ages of globular clusters often exceed estimates of the age of the universe. Age determination of globular clusters 
depends in part on calibration of the absolute magnitude of RR Lyrae stars as a function of period. The 20 brightest RR 
Lyrae stars range from V=7.6 to 10. Parallax measurements are needed at the 1% level, which corresponds to 40 pas for 
RR Lyrae itself. (4) Distances to 90 nearby field subdwarfs cataloged by Carney, 8 ranging in magnitude from 7.2 to 12, 
would calibrate subdwarf luminosities, which are used in fitting the globular cluster main sequences. Parallaxes are 
needed with 30 pas precision. (5) Parallaxes of bright (V = 10) K giants in the galactic disk, accurate to 50 pas. would 
probe the dark matter in our galaxy. 

With 100 pas measurements and a sufficient mission duration, Newcomb could detect Jupiter-sized planets 
around nearby stars. The problem here is that, if the solar system is a "typical" planetary system, signatures large enough 
to be seen will have periods much longer than the 27 month life envisioned for Newcomb. Hipparcos measurements 
might help break the degeneracy between proper motion and a short arc of orbital motion, but this speculation has not yet 
been tested by a sensitivity study. Similarly. Newcomb data could supplement a POINTS mission, especially in the case 
of a planet period long compared to the mission life possibly shortened by equipment failure. 

Gravity Probe B, an experiment which will measure the general relativistic frame dragging due to the spinning 
Earth, needs a bright guide star (Rigel) with proper motion known in an inertial frame to **1 mas/year or better. 

Newcomb could determine the proper motion of Rigel relative to a few bright quasars tied directly to the reference grid. 
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